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The Universe Was Reionized Twice 

Renyue Cen^ 

ABSTRACT 

We show the universe was reionized twice, first at2~15 — 16 and second at 
z ~ 6. Such an outcome appears inevitable, when normahzing to two well deter- 
mined observational measurements, namely, the epoch of the final cosmological 
reionization at z ~ 6 and the density fluctuations at 2; ~ 6, which in turn are 
tightly constrained by Lya forest observations at z ~ 3. These two observations 
most importantly fix the product of star formation efficiency and ionizing photon 
escape fraction from galaxies at high redshift. The only major assumption made 
is that the initial mass function of metal-free. Population III stars is top-heavy. 

To the extent that the relative star formation efficiencies in gaseous minihalos 
with H2 cooling and large halos with atomic cooling at high redshift are unknown, 
the primary source for the first reionization is still uncertain. If star formation 
efficiency in minihalos is at least 10% of that in large halos, then Pop III stars 
in the minihalos may be largely responsible for the first reionization; otherwise, 
the first reionization will be attributable largely to Pop III stars in large halos. 
In the former case, H2 cooling in minihalos is necessarily efficient. We show that 
gas in minihalos can be cooled efficiently by II2 molecules and star formation 
can continue to take place largely unimpeded throughout the first reionization 
period, as long as gas is able to accumulate in them. This comes about thanks 
to two new mechanisms for generating a high X-ray background during the Pop 
III era, put forth here, namely. X-ray emission from the cooling energy of Pop 
III supernova blast waves and from miniquasars powered by Pop III black holes. 
Consequently, H2 formation in the cores of minihalos is significantly induced to 
be able to counteract the destruction by Lyman- Werner photons produced by 
the same Pop III stars. In addition, an important process for producing a large 
number of II2 molecules in relic HII regions of Pop III galaxies, first pointed out 
by Ricotti, Gnedin, & ShuU, is quantified here. It is shown that II2 molecules 
produced by this process may overwhelm the Lyman- Werner photons produced 
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by stars in the same Pop III galaxies. As a result, the Lyman- Werner background 
may never build up in the first place during the Pop III era. 

The long cosmological reionization and reheating history is complex. From 
2; ~ 30 Pop III stars gradually heat up and ionize the intcrgalactic medium, 
completing the first reionization at 2; ~ 15 — 16, followed by a brief period of 
A2; ~ 1, during which the intcrgalactic medium stays completely ionized due 
to sustained ionizing photon emission from concomitant Pop III galaxies. The 
transition from Pop III stars to Pop II stars at 2; ~ 13 suddenly reduces, by 
a factor of ~ 10, ionizing photon emission rate, causing hydrogen to rapidly 
recombine, marking the second cosmological recombination. Prom z ~ 13 to 
z = 6, Compton cooling by the cosmic microwave background and photoheating 
by the stars self-regulate the Jeans mass and the star formation rate, giving rise 
to a mean temperature of the intergalactic medium maintained nearly at constant 
of ~ lO'' K. Meanwhile, recombination and photoionization balance one another 
such that the intergalactic medium stays largely ionized during this stage with 
nuii/nu > 0.6. Most of the star formation in this period occurs in large halos 
with dominant atomic line cooling. 

We discuss a wide range of implications and possible tests for this new reion- 
ization picture. In particular, the Thomson scattering optical depth is increased 
to 0.10 ± 0.03, compared to 0.027 for the case of only one rapid reionization at 
z = 6. Upcoming Microwave Anisotropy Probe observation of the polarization 
of the cosmic microwave background should be able to distinguish between these 
two scenarios. 

Subject headings: cosmology: theory — intergalactic medium — interstellar 
medium — supernova — reionization 

1. Introduction 

The conventional view is that the universe becomes reionized at some point in the red- 
shift range 2; = 6 — 10, when the UV emission rate in galaxies with virial temperatures greater 
than ~ lO^K (where hydrogen atoms are efficient coolants) exceeds the overall recombination 
rate (e.g., Barkana & Loeb 2001; Madau 2002). At the time of this writing the redshift range 

has been narrowed to a point at ^ ~ 6, as suggested by recent observations of high redshift 
quasars from the Sloan Digital Sky Survey (SDSS) (e.g.. Fan et al. 2001; Becker et al. 2001, 
Barkana 2001; Cen & McDonald 2002; Litz et al. 2002). 

In this paper we present a new scenario. In the context of the standard cold dark 
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matter cosmological model we show that the universe was first reionized at z = 15 — 16 by 
Population III (Pop III) stars and second reionized at ^ = 6. Following the first reionization, 
the transition from Pop III stars to Pop II stars occurs. At this point photoionization becomes 
insufficient to counterbalance the rapid recombination process and the intergalactic medium 
(IGM) recombines to become opaque to Lja and ionizing photons, again! As time progresses, 
with increased density fluctuations and the nonlinear mass scale star formation rate gradually 
picks up. At z ~ 6 the global star formation rate exceeds the global recombination rate 
and the universe is completely reionized for the second time, and for all. Prom the first 
cosmological reionization through the second cosmological reionization the mean temperature 
of the IGM is maintained at ~ 10^ K balanced between Compton cooling and photo-heating 
and hydrogen is more than half ionized balanced between recombination and photoionization. 

We show that this new reionization picture is inevitable, as long as Pop III initial mass 
function (IMF) is top-heavy, thanks to a new, powerful constraint placed on the product of 
star formation efficiency and ionizing photon escape fraction from galaxies at high redshift, 
which otherwise would be unknown. This constraint comes from two solid pieces of obser- 
vations: 1) the universe is required to be reionized at ^ ~ 6 and 2) the density fluctuation 
in the universe at 2; ~ 6 is well determined by the same small-scale power traced by the 
observed Lja forest observed at 2; ~ 3. 

While it is clear that the ionizing sources for the second cosmological reionization are 
stars in large galaxies with efficient atomic line cooling, the ionizing sources for the first 
cosmological reionization are uncertain to the extent that we do not know what the relative 
star formation efficiencies in large halos and minihalos are. We define "minihalos" as those 
whose virial temperature is less than ~ 8 x lO^K where only H2 cooling is possible in 
the absence of metals, and large halos as those with virial temperature above ~ 8 x lO^K 
capable of cooling via atomic lines. We show that, if star formation efficiency in minihalos 
with H2 cooling is not more than a factor of ten less efficient than in large halos with atomic 
cooling, then Pop III stars in minihalos may be largely responsible for the first reionization. 
Conversely, the first reionization may be attributable largely to Pop III stars in large halos. 

In order to enable efficient H2 cooling in minihalos it is necessary to maintain an ade- 
quate level of H2 molecule fraction within their cores. As is well known, H2 molecules are 
fragile and easily destroyed by photons in the Lyman- Werner (LW) bands (11.18 — 13.6eV; 
Field et al. 1966; Stecher & Williams 1967), to which the universe is largely transparent. 
Primeval H2 molecules have been completely destroyed well before enough ionizing photons 
are produced to reionize the universe (Gnedin & Ostriker 1997; Haiman, Rees, & Loeb 1997; 
Tegmark et al. 1997). Therefore, an adequate production rate of H2 molecules is required to 
counteract the destruction rate to keep H2 fraction at a useful level. A sufficiently high X- 
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ray background at high redshift could serve as a requisite catalyst for forming H2 molecules 
by deeply penetrating into and generating a sufficient number of free electrons in the cores 
of minihalos (Haiman, Rees, & Loeb 1996; Haiman, Abel, & Rees 2000; Ricotti, Gnedin, & 
ShuU 2001, RGS hereafter; Glover & Bland 2002). 

We put forth two new mechanisms for generating a high X-ray background during the 
Pop III era. We point out that at ^ = 13 — 20 Pop III supernova remnants and miniquasars 
powered by Pop III black holes are efficient X-ray emitters. The much higher density of the 
interstellar medium at high redshift results in a much more rapid and earlier cooling phase 
of the supernova blast waves, emitting a large fraction of the cooling energy in X-ray energy 
(~ IkeV). [In contrast, supernova blast waves in the local galaxies are known not to cool 
efficiently until their temperatures have dropped well below the X-ray regime (e.g.. Cox 1972; 
Chevaher 1974,1982)]. Miniquasars powered by Pop III black holes of mass ~ 10 — 100 Mq 
are expected to emit a large fraction of the radiation in X-rays at > 1 keV. The combined X- 
ray emission from these two sources is sufficient to enable efficient H2 formation and cooling 
in the cores of minihalo galaxies in the redshift period in question. 

In addition, an important process for producing a large number of H2 molecules in relic 
HII regions of Pop III galaxies is quantified. It is shown that H2 molecules produced by this 
process may overwhelm the Lyman- Werner photons produced by stars in the same Pop III 
galaxies. As a result, the Lyman- Werner background may have never built up in the first 
place during the Pop III era. In combination, we suggest that star formation in gas-rich 
minihalos can continue to take place largely unimpeded throughout the first reionization 
period. 

The outline of this paper is as follows. Two new sources of X-ray radiation, namely, 
the Pop III supernova cooling radiation and X-ray emission from miniquasars are quantified 
in §2. In §3 we quantify the II2 formation in relic HII regions of Pop III galaxies. In §4 
we compute in detail the evolution of the IGM from z ^ 20 to z = 6. We discuss possible 
implications and significant consequences and tests of this new scenario in §5, and conclude 
in §6. Throughout a spatially flat cold dark matter cosmological model with Qm — 0.25, 
Qb = 0.04, A = 0.75, Ho = 72km/s/Mpc and erg = 0.8 is adopted (e.g., BahcaU et al. 2002). 

2. X-ray Emission Prom Pop III Supernovae and Miniquasars, and Moleculcir 

Hydrogen Formation in Minihalos 

The important role that H2 molecules play in the collapse of gas clouds to form Pop 
III stars and the initial mass function of Pop III stars have a long history of investigation 
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(Saslaw & Zipoy 1967; Hirasawa, Aizu, & Taketani 1969; Takeda, Sato, & Matsuda 1969; 
Hutchins 1976; Silk 1977; Hartquist & Cameron 1977; Shchekinov & Edelman 1978; Yoshii 
& Sabano 1979; Tohline 1980; Carlberg 1981; Lcpp & Shull 1984; Yoshii & Saio 1986; Lahav 
1986; Stabler 1986; Shapiro & Kang 1987; Uchara et al. 1996; Haiman, Thoul, & Loeb 1996; 
Padoan, Raul, & Jones 1997; Nakamura & Umemura 1999,2001,2002; Larson 1995,2000; Abel 
et al. 1998; Abel, Bryan, & Norman 1999,2000,2002; Bromm, Coppi, & Larson 1999, 2002; 
Fuller & Couchman 2000; Machacek, Bryan, & Abel 2001). Many authors have examined 
the possibility of Pop III stars in minihalos reionizing the universe (e.g., Couchman & Rees 
1986; Fukugita & Kawasaki 1994; Haiman & Loeb 1997) or partially reionizing the universe 
(e.g., Gnedin & Ostriker 1997; Gnedin 2000a; RGS). The primary difficulty for Pop III stars 
in minihalos to reionize the universe has by now been realized to be the destruction of H2 
molecules by photons in the LW bands produced by the same Pop III stars, long before the 
completion of reionization process (Gnedin & Ostriker 1997; Haiman, Rees, & Loeb 1997; 
Tegmark et al. 1997); H2 photodissociation time becomes shorter than the Hubble time 
when the ionizing radiation intensity at the Lyman limit reaches ~ 10^^^ erg/cm^ /sec/hz/sr 
at the redshift of interest, about three orders of magnitude below what is required to ionize 
the universe. 

With regard to the initial mass function of Pop HI stars a new picture is emerging 
from a number of recent theoretical studies of the collapse of primordial gas clouds at high 
redshifts, induced by H2 coohng. At redshift 2; ~ 10 — 30 primordial gas clouds are shown 
to collapse to form very massive stars with mass M > 100 Mq (Abel et al. 2002; Bromm 
et al. 2002; Nakamura & Umemura 2002). This outcome is basically determined by the 
Jeans mass of the collapsing cloud, involving a complicated interplay between cooling and 
fragmentation. For Pop HI galaxies at the redshift in question (2; ~ 13 — 20) the initial (i.e., 
interstellar) gas density is less than the threshold value of 10^ cm^"^ identified by Nakamura 
& Umemura (2001,2002) and fragmentation is to occur at low density with fragment mass 
of ~ 100 Mq. We adopt this new theory for Pop HI stars, called Very Massive Stars (VMS). 
However, we note that the main conclusions drawn in this paper will remain unchanged as 
long as the initial mass function of Pop HI stars is substantially top-heavy (e.g., Umeda & 
Nomoto 2003). 

We point out two new important sources for X-ray emission in this section, capable of 
creating a high X-ray background at high redshift to enable the formation of enough H2 
molecules to induce continuous cooling and star formation in minihalos during the Pop HI 
era. 
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2.1. X-ray Emission from Pop III Supernova Remnants 

We use a simple analytic means to estimate the X-ray emission from Pop III supernova 
remnants assuming that a spherical supernova blast wave propagates into a uniform-density 
interstellar medium. A spherically expanding supernova remnant at sufficiently late times 
follows the simple, adiabatic, self-similar Sedov-Taylor solution (Sedov 1959; Taylor 1950; 
Shklovsky 1968; Cox 1972; Ostriker & Cowie 1981). Specifically, the shock radius (Rs), 
shock velocity (Vg) and postshock temperature (Tg) (for 7 = 5/3) obey: 

R, = 21.7tT{E,2/ny/^ pc, (1) 

= 83Mf^\E52/nY^'' km/s, (2) 

T, = 8.45 X 10%'^\E,2/nf^' K, (3) 

where E^2 is the explosion energy in 10^^ ergs, is the time elapsed since the onset of 
explosion in 10"^ yr and n is the density of the interstellar medium in cm~^. 

The Sedov-Taylor phase ends with a rapid cooling phase, resulting in the formation of a 
thin dense shell (e.g., Ostriker & McKee 1988). This radiative cooling phase of the shocked 
gas sets in abruptly approximately when the cooling time is equal to the time elapsed, i.e., 
when 

where 100% hydrogen is assumed in computing the internal energy and the electron density 
Ue for simplicity (but not for the cooling function A; see below); k is the Boltzmann's con- 
stant; He = 4:12 is the postshock electron density; A(Te) is the volume cooling function for a 
primordial plasma with 76% hydrogen and 24% in helium by mass, taken from Sutherland 
& Dopita (1993). Since the blast waves of the first generation of supernovae propagate into 
a primordial gas free of metals, metal cooling is non-existent. More importantly, complex 
processes due to dust need not to be considered; in particular. X-ray emission would not 
suffer significantly from dust absorption (Ostriker & Silk 1973; Burke & Silk 1974; Draine 
& Salpeter 1979; ShuU 1980; Wheeler, Mazurek, Sivaramakrishnan 1980; Tielens et al. 1987; 
Draine & McKee 1993). Combining equations (3,4) and solving for Ts as a function of n give 
results shown in Figure 1, for two cases with (£^52, Z/ Zq) = (5.0, 0.0) and (0.1, 1.0), respec- 
tively, where the former may be appropriate for supernovae/hypernovae resulting from very 
massive Pop III stars (Woosley & Weaver 1982; Ober et al. 1983; Bond et al. 1984; Nakamura 
et al. 2001; Heger & Woosley 2002) and the latter for present-day normal supernovae. 

A critical point to note is that the density of the general interstellar medium should scale 
with [1 + zY, being significantly higher at high redshift than that of the local interstellar 
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Fig. 1. — The solid curve shows the temperature of the coohng shock front for a super- 
nova/hypernova from a very massive Pop III star with initial energy of 5 x 10^^ erg as a 
function of the interstellar (metal- free) gas density (bottom x-axis) or I + 2; (top x-axis). The 
cooling function assumes a gas of primordial composition. The vertical tick hanging from 
the top X-axis indicates redshift z = 15.5 and the vertical tick standing on the bottom x-axis 
indicates redshift z = 0. Cooling functions are taken from Sutherland & Dopita (1993). As 
a reference, the dashed curve is for a gas with solar metallicity and a supernova explosion 
energy of 10^^ erg, more appropriate for a normal supernova. 
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medium. This assumption should hold in the context of cosmological hierarchical structure 
formation for the following reasons. First, the mean gas density scales with (l + z)^. Second, 
halos at low and high redshift in cosmological simulations show similar properties when 
density and length arc measured in their respective comoving units (e.g., Navarro, Frenk, 
& White 1997; Del Popolo 2001). Third, the spin parameters (i.e., angular momentum 
distribution) of both high and low redshift halos have very similar distributions peaking at 
a nearly identical value A ~ 0.05 (Peebles 1969; White 1984; Barnes & Efstathiou 1987; 
Ueda et al. 1994; Steinmetz & Bartelmann 1995; Cole & Lacey 1996; Bullock et al. 2001; 
Cen et al. 2003). Thus, coohng gas in galaxies at low and high redshift should collapse 
by a similar factor before the structure becomes dynamically stable (e.g., rotation support 
sets in), resulting in interstellar densities scaling as (1 + z)^. Direct simulations (Abel et al. 
2000,2002; Bromm et al. 2002) suggest a gas density of 10"^ — lO'^cm^'^ by the end of the initial 
free fall for minihalos at 2; ~ 20, verifying this simple analysis. The same argument was given 
in Kaiser (1991) on a somewhat larger scale for clusters of galaxies. Peebles (1997), using 
the same scaling of density with respect to redshift, suggested that an early epoch of galaxy 
formation may be favored. Observationally, the higher density in bulges of spiral galaxies 
or elliptical galaxies are consistent with their earlier formation epochs. In Figure 1 we have 
used n{z) — no{l + z)^, with no — lcm~^ being the density of the local interstellar medium, 
to translate the bottom x-axis (n) to the top x-axis {1 + z). 

In Figure 1 we see that, ai z — 15.5 (or n — 4492cm~^), the blast wave enters the rapid 
coohng phase at a postshock temperature shghtly above 1.0 x lO^K. Most of the energy is 
radiated away during the brief coohng phase (Falle 1975,1981) with photon energy hu ~ 
kTs. [Subsequent evolution of the cooling shell will be subject to various instabilities (e.g.. 
Chevalier & Imamura 1982; Vishniac 1983; Bcrtschingcr 1986; Ciofii, McKcc & Bertschinger 
1988). Following that, the evolution enters snowplow phase driven by the pressure of a still 
hot interior gas (McKcc & Ostrikcr 1977; Ostrikcr & McKcc 1988).] For a primordial gas 
cooling at 1.0 x lO'^K, it is found that (40%, 31%, 18%) of the instantaneously radiated energy 
is at photon energies above (0.8, 1.0, 1.5)keV, respectively. Clearly, a significant amount of 
the total energy of the supernova blast wave will be turned into X-ray photons. As we will 
see later, 2; = 15 — 16 will be identified as the redshift of the first cosmological reionization by 
Pop III stars. At higher redshift, the emitted photons from supernova remnant shell cooling 
would be still harder. 

As a consistency check we find that, for n — 4492cm~^, at the onset of the rapid cooling 
phase, the elapsed time is t^ad — 857 yr, the shell radius is Trad — 2.09pc, the swept-up 
interstellar medium mass is Mrad = 4169 M©. The fact that Mrad is much larger than the 
mass of the supernova ejecta (Mej ~ 100 M0) guarantees the Sedov- Taylor solution for the 
regime in question (until cooling sets in). The fact that the shell radius Trad is much smaller 
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than the size of galaxies (> lOpc) and Mrad is much smaller than the total baryonic mass 
in the galaxies (> 10^ M©) in question indicates that the blast wave at the cooling time 
is still sweeping through high density interstellar medium, as was assumed. As another 
consistency check, the time required for electrons and ions to reach temperature equilibrium 
is teq = 1.4 X 10'^£'52^'^n^^/'^ = 162 yr for the case considered, which is much shorter than 
trad, indicating that electrons can radiate away the shock heated thermal energy and that 
Sedov-Taylor phase is vahd. In contrast, in the star-burst model for AGN (Terlevich et 
al. 1992), the proposed supernova remnants interact with a much higher density medium 
{n ~ lO^cm"^) such that the Sedov-Taylor phase is never reached due to extremely rapid 
cooling before thermalization. 

In brief, a large fraction of supernova explosion energy in Pop III galaxies at 2; ~ 13 — 20, 
possibly as large as /x ~ 0.30, is shown to be converted into X-ray photons with energy 
greater than IkeV. The X-ray background produced by this process will be shown below to 
be able to play a positive feedback role in the formation of and cooling by H2 molecules in 
minihalos. 

Let us now proceed to identify the characteristics of the composite spectrum of the 
background radiation field produced by both emissions from the VMS and the thermal 
emission from VMS supernova blast waves. A relevant parameter for our purpose is the 
ratio of the energy in the LW bands {hu = 11.18 — 13.6eV) to the energy in photons with 
hu > 1 keV, \E'in = £'Lw/-E'>iA:eV'- Photons in the LW bands (11.18 — 13.6eV) are primarily 
responsible for photodissociating H2 molecules and are not heavily absorbed by intervening 
intergalactic atomic hydrogen (except for the saw-tooth modulation by the atomic Lyman 
hue series; Haiman, Abel, & Rees 2000; HAR hereafter). Hard X-ray photons at hu > 
1 keV, on the other hand, are also largely unabsorbed by intervening atomic hydrogen 
and helium and capable of penetrating deeply into the cores of minihalos to produce free 
electrons through both direct photo-ionization and secondary photo-electron ionization. The 
abundance of H2 molecules in the cores of minihalos is primarily a result of the competition 
between the two. In contrast, photons with energy in the range 13.6eV— IkeV are heavily 
absorbed by atomic hydrogen and helium prior to complete reionization of the universe and 
thus have little effect on the formation of H2 molecules in the cores of minihalos (byt see 
RGS for a positive feedback due to ionizing photons on the surface of the expanding HII 
regions and in relic HII regions). We may write approximately as 

*m ~ ^-7 , (5) 

where /lw is the fraction of energy in the LW bands emitted by VMS; c is the speed of light; 
-^VMS is the characteristic mass of VMS (~ 200 M0; see §3); E^^ is the explosion energy of a 
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typical VMS; ^imf < 1 is inserted to take into account the effect that some VMS (for example, 
non-rotating stars with M > 260 Mq; Rakavy, Shaviv, & Zinamon 1967; Bond, Arnett, & 
Carr 1984; Glatzel, Fricke, & El Eid 1985; Woosley 1986) do not produce supernovae (but 
collapse wholly to black holes). The supernova explosion energy, E^x, released by VMS in 
the mass range 140 — 260 M© is approximately 10^^ — 10^^ erg (Woosley & Weaver 1982; Ober 
et al. 1983; Bond et al. 1984; Nakamura et al. 2001; Heger & Woosley 2002). The VMS has 
approximately a blackbody radiation spectrum with an effective temperature of ~ 10^'^ K 
(Tumlinson & ShuU 2000; Bromm et al. 2001) for which it is found that /lw = 1-5 x 10-^, 
yielding 

*n,=2.5«-(A)-'(^^J^)-'(^). (6) 

Is it small enough to enhance H2 formation in minihalos? Let us compare to the results of 
a systematic study of the effect of X-ray photons on the cooling of minihalo gas by HAR. 
HAR investigated the effect using a. oc background radiation spectrum, with an upper 
cutoff at 10 keV. Using ^ to parameterize their results, the findings of HAR are that, when 

^ 1d(13.6/11.18) 

*™ ^ ..1.(10/1) = 

(for Ex — 0.016), H2 formation is enhanced and cooling by H2 is sufficient to allow for the 
gas in the cores of minihalos with virial temperature Ty > lOOOK to collapse to form stars. 
We note that in evaluating \E'har above we have taken into account a factor of ~ 6 underesti- 
mation of ^HAR in the original work of HAR due to a factor of 6 underestimation of the LW 
photon absorption cross section by H2 molecules (Haiman 2002, private communications; 
RGS; Glover & Brand 2002). The correction factor 6 is a lower hmit in the optically thin 
case for H2 molecules for the cores of minihalos. In the optically thick limit, a substantially 
larger correction factor needs to be applied (i.e., ^'har would be substantially larger than 
5.1), although the exact effect would require a re-calculation of the results in HAR. It thus 
appears that the X-ray emission from Pop HI supernova remnants alone may be able to 
produce enough X-ray photons relative to the number of photons in the LW bands such 
that production rate of H2 molecules dominates over the destruction rate in most of the 
minihalos. 



2.2. X-ray Emission from Miniquasars Powered by Pop III Black Holes 

Another direct consequence of Pop HI VMS is an inevitable production of a significant 
amount of black holes. Pop HI VMS more massive than 260 M0 would eventually implode 
carrying the entire mass to form black holes (Rakavy, Shaviv, & Zinamon 1967; Bond, Arnett, 
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fe Carr 1984; Glatzel, Pricke, & El Eid 1985; Woosley 1986) without producing explosions. 
Pop III VMS less massive than 140 M0, after about three million years of luminous life, 
would explode as supernovae or hypernovae leaving behind black holes of mass ~ 10 — 50 Mq 
(Hcger & Woosley 2002). These Pop 111 black holes could accrete gas from surroundings to 
shine as miniquasars. The likelihood of gas accreting onto these black holes is at least as 
high as in their lower redshift counterparts since the gas density is higher and halos tend to 
have somewhat lower spins at high redshift. In the case of the very massive Pop III black 
holes (M > 260 Mq) formed without explosion, gas may be ready to accrete immediately, 
as the surrounding gas has not been blown away. The dynamics of the black holes in these 
environments in the context of cosmological structure formation is a complex subject and 
can only be treated in separate works. Here we make the assumption that these black holes 
would accrete gas and grow. 



The characteristic gas temperature of a disk powered by accretion onto a black hole is 
(Rees 1984): 



where M2 is the black mass in IO^Mq. For Mbh = (2OO,2O)M0, we have Te = (1.3 x 
10^, 2.4 X 10^) K. While the spectral energy distribution (SED) of quasars powered by su- 
permassive black holes is known to contain a significant fraction in X-rays (e.g., Elvis et al. 
1994), the miniquasars at high redshift powered by much smaller black holes will be con- 
spicuous in X-rays, probably emitting predominantly in X-ray band from both thermal and 
nonthermal emission. A somewhat more quantitative argument may be made as follows. The 
SED of observed quasars powered by supermassive black holes of mass ~ 10^ Mq contains a 
substantial amount of energy in the X-rays but the largest concentration of energy appears 
to peak at ~ 12 — 13eV (UV bump; Elvis et al. 1994) barring the unknown gap between 
UV and X-ray. Under the reasonable assumption that the SED peak frequency scales with 
the characteristic temperature Te, then the peak frequency for Pop III black hole powered 
miniquasars would be shifted (by a factor of ~ 100 compared to quasars) to hi/ ~ IkeV. 

Let us now compute the X-ray background produced by miniquasars powered by Pop III 
black holes. We put this quantitatively in the context of the ratio of energy in LW photons to 
that of X-ray photons. First, an order-of-magnitude estimate. If a miniquasar has accreted 
a fraction, rjQ, of the initial mass of the black hole, the radiated energy will be 



Te = 1.3 X lO^Ma" ' K, 



(8) 



Eq = riqac^MBH, 



(9) 



where a is the radiative efficiency. Then, the ratio ^ may be written as 



0.007cVm/LW 



(10) 



O.lao.iC^Tyg/in/BH/xray/c 



esq 



- 12 - 



where /m is the fraction of mass collapsed into Pop III stars; /bh is the fraction of mass in 
fill that ends in black holes; ao.i is radiative efficiency of the miniquasars in units of 0.1; /xray 
is the fraction of energy radiated by miniquasars with photon energy greater than 1 keV; 
/esq is the X-ray escape fraction from miniquasars; we have used /lw = 1-5 x 10^^. The 
value of ao.i is thought to be close to unity (Rees 1984) and the value of /xray is expected to 
be of order unity. The fraction of mass in Pop 111 stars ultimately collapsing to black holes 
/bh is also of order unity, if either Pop III stars are very massive ~ 100 Mq or » 1 Mq. The 
accretion fraction rjQ would depend on the gas accretion rate and the formation time of the 
black hole; parameterization by Equation (10) would be fairly accurate if black holes formed 
long time ago. But a more accurate way to compute the X-ray energy from miniquasars is 
to integrate over time over all miniquasars by parameterizing each radiating at ^Edd times 
the Eddington luminosity. The results are shown in Figure 2, where the solid and dotted 
curves show with radiative efficiency of a = (0.1,0.03), respectively. The global Pop 111 
star formation rate is computed using Press-Schcchtcr (1974) formula considering all halos 
with virial temperature greater than 1000 K. The reason that the dotted curve with lower 
radiative efficiency (a = 0.03) lie below (i.e., emits more X-ray photons than) the solid curve 
is that the black hole mass in this case grows at a faster rate, being inversely proportional 
to the radiative efficiency. 

Another useful way to estimate the growth of mass in black holes is to follow the local 
observed black hole mass-bulge mass relation, Mbh ~ 0.006Mbuige, (Magorrian et al. 1998). 
If the same ratio of the mass in black holes to bulge mass, which at high redshift may be 
equated to the total stellar mass, holds at high redshift, then we can readily compute the 
energy hberated by gas accretion onto Pop III black holes as 

vn _ 0.007cViii/lw _ 1 //esq\ ._i ._i 

0.006 X O.lao.lC^/m/BH/xray/esq V 0-2 / ^ 

The horizontal dashed line in Figure 2 indicates computed this way with radiative effi- 
ciency of q; = 0.1 and X-ray escape fraction /esq = 0.2. 

From Figure 2 it is evident that ^'q is significantly smaller than "ifnAR = 5.1 (Equation 
7), at z < 20, required to induce sufficient H2 coohng in minihalos (HAR). The higher value 
of at higher redshift {z > 20) is due to the fact that the age of the universe becomes 
much shorter than the Eddington time of ~ 4 x 10^ yr and the black holes have not had 
enough time to accrete a substantial amount of gas; the age of the universe is 1.9 x 10^ yr 
at ^ = 20 (for H ^ 70 km/s/Mpc, Qm = 0.25, A = 0.75). 

In summary, two distinctive mechanisms, namely supernova remnants and miniquasars, 
each appears to be able to generate enough X-ray emission at 2; ~ 13 — 20 to provide positive 
feedback on star formation in subsequent (other) minihalos by making a sufficient number 
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Fig. 2. — The solid and dotted curves show \1/q as a function of redshift, by integrating over 
time the miniquasar luminosity radiating at Eddington luminosity with radiative efficiency 
of q; = (0.1,0.03), respectively. The horizontal dashed line indicates the case if one adopts 
the local observed black hole mass-bulge mass relation (Magorrian et al. 1997) with radiative 
efficiency of a = 0.1. The horizontal long-dashed fine indicates the required value (HAR), 
below which H2 cooling is efficient, would scale with ^Edd and f^nfxLy^ both of 
which we have assumed value of unity. We also assume that 20% (/esq = 0.2) of the soft 
X-ray radiation (~ 1 keV), averaged over all miniquasars, escapes to the intergalactic space; 
this value is based on observations of AGN and quasars at lower redshift (Brandt 2002, 
private communications). We have ignored cosmological redshift and dilution effects in this 
illustration. 



-14- 



of X-ray photons relative to the number of destructive photons produced in the LW bands. 
The combination of the two should ensure that enough X-ray radiation is produced as a 
result of Pop III star formation. In addition, there may be other, significant X-ray emission 
mechanisms such as inverse Compton emission (e.g., Hogan & Layzer 1979; Oh 2001) or 
massive X-ray binaries (e.g.. Bookbinder et al. 1980; Helfand & Moran 2001). Moreover, 
additional positive feedback mechanisms, such as that proposed by Ferrara (1998) due to 
enhanced H2 formation in the supernova cooling shells, that put forth by RGS from the 
enhanced H2 formation at the surfaces of Stromgren spheres of individual Pop III galaxies 
and that quantified in the next section (§3) due to H2 formation in relic H II regions produced 
by Pop III galaxies will further help promote H2 formation. Taking all the processes together, 
it appears highly likely that the chief obstacle to continuous H2 formation and cooling in 
minihalos is removed. 

Since the assumption that X-ray emission produced by these two processes related to 
Pop III star formation can produce positive feedback to subsequent star formation is impor- 
tant, it is warranted to have independent checks. We will compare our results to Glover & 
Brand (2002; GB hereafter). We compare to their X-ray emission model due to inverse Comp- 
ton, which has luminosity density Lx = 7.7 x 10^^(z/o/^)/c erg s^-'^Hz^^( MQyr^-*^)^-*^, where 
Uq = IkeV and /e is the fraction of supernova energy transferred to rclativistic electrons. 
Integrating Lx over a range 1 — lOkeV we obtain Ex — 4.28 x 10^^ /e erg s"^ ( MQyT^^)^^ . The 
luminosity at the LW bands in GB is Llw = 1-1 x lO^^erg s~^ Hz~^( M0yr~^)~^, which after 
integration over LW bands gives E'lw — 6.4 x lO^^erg s"^ (Mgyr"^)"^. Taking the ratio of 
the above two luminosities yields ^^gb = ^'LwZ-E'x = 150, for the most optimistic model that 
they use with /e = 0.1. For this model GB find that the critical virial temperature of halos 
where cooling is sufficient is lowered from 8 x lO'^K to about (2 — 3) x lO^K at 2; ~ 17 (see 
Figure 6 in GB), a significant effect. While even this model is already sufficient to enable 
molecular cooling in halos with T„ > 2000K to reionizc the universe at 2; ~ 17 (sec Figure 
5 below), it is expected that the much smaller ^'m computed above for Pop III stars in our 
model would further drive the critical virial temperature to much lower values and enable 
star formation in almost all minihalos, consistent with the conclusions reached by HAR. 

3. Formation of H2 Molecules in Relic HII Regions Produced by Pop III 

GalaLxies 

While the positive feedback due to a high X-ray background produced by Pop III super- 
novae and miniquasars powered by Pop III black holes seems sufficient to sustain continuous 
gas cooling in gaseous minihalos, we quantify another important source for production of 
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H2, which was first pointed out by RGS. 

Each Pop III galaxy creates an HII region (Stromgren sphere) around it. Since the 
hfetime of a Pop III galaxy of 3 x 10^ yr is comparable to and somewhat shorter than the 
hydrogen recombination time, tree = 1.3 x lO^yr at 2; ~ 20, we may assume for simplicity, 
without causing substantial errors, that the size of the Stromgren sphere is just determined 
by the number of ionizing photons emitted, A^ion- We follow the evolution of such a relic HII 
region after the death of a Pop III galaxy, in the absence of any external radiation field. We 
solve a set of rate equations and energy equation for relevant species in an H II region free 
of metals. The relevant reaction coefficients are taken from Abel et al. (1997). 

Figure 3 shows the fraction of II2 molecules as a function of redshift for four H II regions 
formed aX z = (17, 18, 19,20), evolved up to redshift z = 15, which will be identified as the 
first rcionization epoch by Pop III stars (§4). We see that the H2 fraction in H II regions 
shows a sharp rise at an early time and flattens out at later times. The reason for this 
asymtotic convergence may be understood by examining Figure 4, where the temperature 
evolution of the gas in the region is shown. H2 formation rate decreases with decreasing 
temperature (Abel et al. 1997), causing the total H2 abundance to flatten out below 1000 K. 
On the high temperature end, since H2 formation rate only becomes signiflcant after the 
temperature has dropped below about 7000 K, we see a slight delay {Az ~ 1.5) in the ascent 
of the H2 abundance in Figure 3. 

Let us now assess the signiflcance of the H2 molecules formed in these rehc H II regions. 
A critical ratio here is the number of H2 molecules produced over the number of LW photons 
produced, denoted as -Rreiic- For A^ion ionizing photons emitted by Pop III stars, there will 
be approximately A^hii = A^ion hydrogen atoms ionized in the H II region. At the same time, 
the number of LW photons produced by the same stars is 1.5 x lO^^A^ion- From Figure 3, 
we take an approximate average nH,^/nH = 2 x 10^^ for the abundance of H2 molecules 
produced hj z = 15. Then, we obtain -Rreiic = 0.13. Since about 10 LW photons are needed 
to photodissociate one H2 molecule (RGS; Glover & Bland 2002), the break-even point for 
the ratio of the number of H2 molecules produced over the number of LW photons produced 
is about 0.1, i.e., no net molecules would be produced. It thus appears that H2 molecules 
produced in the relic H II regions produced by Pop III stars may be sufficient to offset 
the production of LW photons by the same stars. Consequently, LW radiation background 
may never be able to build up during the Pop III era. However, a more careful calculation, 
taking into account the density fiuctuations among others, is needed to precisely quantify this 
positive feedback effect. But we defer such a calculation to a future work. Nevertheless, it 
seems possible that this positive feedback in relic H II regions may be more effective than that 
proposed by Ferrara (1998) due to enhanced H2 formation in the supernova cooling shells 




Fig. 3. — shows the fraction of H2 molecules as a function of redshift for four relic HII 
regions formed at z = (17,18,19,20), respectively. We assume a reasonable and perhaps 
conservative clumping factor of 15 for HII regions in the vicinity of Pop III galaxies. 
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Fig. 4. — shows the evolution of temperature for four H II regions formed a.t z = 
(17, 18, 19,20), respectively. 
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and that put forth by RGS from the enhanced H2 formation at the surfaces of Stromgren 
spheres of individual Pop III galaxies. 

In summary, it appears that LW radiation background may have never been able to 
build up to destroy H2 molecules in subsequent minihalos. In the event that LW radiation 
background was able to build up, its intensity would have been substantially reduced. The 

fact that the X-ray background produced by Pop III supernovae and miniquasars powered by 
Pop III black holes have already been capable of countering an undiminished LW radiation 
background (i.e., by ignoring all the H2 molecules produced by these positive feedback 
processes) leads us to conclude that subsequent star formation in minihalos will not be 
hindered by previous star formation in the Pop III era. 



4. The Process of Cosmological Reionization 

4.1. Computational Method 

We will examine the evolutionary history of the IGM from a very early redshift {z » 20) 
until the universe is completely reionized at 2; ~ 6, as suggested by recent observations. 
Rather than taking a brute-force approach, which will be deferred to a later work, we will 
use a new, fast Monte Carlo-hke approach to economically sample parameter space. This 
method, we hope, should be capable of capturing the essential physical processes involved. 

At any time the gas in the universe consists of N distinct regions in the two-dimensional 
phase space specified by {Ti,yi), where Tj and i/i are the temperature and neutral hydrogen 
fraction, respectively, and N is a varying number. Each phase space region i contains 
mass fraction of the total gas in the universe and its mean gas density is assumed to be 
equal to the global mean; the sum of all /j is unity. For each region i we solve a combined 
set of equations simultaneously to follow the coupled evolution of gas and star formation: 

dfi dNion 1 



dt dt ny 

dTi Yii(Ti, yi, Ci, z 



fi, (12) 

(13) 



dt 3k 

^ = aQri(l - yif - (3Qnyi(l - y,), (14) 

where n is the global mean hydrogen number density; y is the global average of the neutral 
hydrogen fraction; Ej is the net cooling rate per unit mass (including Compton cooling, 
atomic line cooling, recombination cooling, photoheating due to ionization and adiabatic 
cooling due to cosmic expansion), which is a function of T^, y^, clumping factor Ci (see 
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below) and redshift z] k is the Boltzmann's constant; a is the case B hydrogen recombination 
coefficient; /3 is the hydrogen coUisional ionization coefficient; dNif^/dt is the number of 
ionizing photons emitted per unit volume per unit time, equal to 

dNian _ c*euvc^pfes dip 
dt ~ 13.6eV dt ^ ^ 

where c* is the star formation efficiency; ejjv is the fraction of energy (in units of rest mass of 
stars formed) turned into hydrogen ionizing photons; p is the global mean gas mass density; 
c is the speed of light; fes is the ionizing photon escape fraction from galaxies; ^ is the 
global rate of fraction of gas formed into stars, equal to 

dip ^rdlpi . s 

1=1 

where ^ is the rate of the fraction of gas formed into stars in region i. At each time step, 
we use the Press-Schechter (1974) formula: 

V^,(>M,^)=erfc( f ) (17) 

to compute the fraction of gas, ipi, that is in halos in region i; where erfc() is the complimen- 
tary error function and cr(M, z) is the rms density fluctuation smoothed over a region of mass 
M at z; 5c is a constant equal to 1.69 indicating the amplitude of the density fluctuations 
on a top-hat sphere that collapses at z (Gunn & Gott 1972). The threshold halo mass M 
is basically the Jeans mass, given the IGM temperature of region i; i.e., only those halos 
whose mass is greater than the Jeans mass at that time will accrete gas and contribute to 
the star formation at that time step. In practice, however, the threshold mass M depends 
on the history of the gas involved as well as non-trivial effects due to dynamically dominant 
dark matter and is thus best determined by detailed simulations. We will use the empirically 
determined "filter mass" formula from full hydrodynamic simulations by Gnedin (2000b): 

Mf = 1.0 X 10««:(-^)-V2 ( J^) ( l±f ^ , (18) 

where at the relevant redshift range, 2; ~ 6 — 20, the empirically determined constant k is 
found to be 0.5 — 1.6. The results obtained in §4.2,4.3 do not depend on k, sensitively in the 
indicated range of k. 

As has been vahdated in §2,3 at the Pop III era, star formation occurs in all halos where 
gas is able to collect, not limited to large halos where atomic hydrogen coohng becomes 
efficient. It is, however, unclear whether H2 cooling would still be efficient in the Pop II era. 
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We show that even in the absence of the X-ray background, a high residual ionization fraction 
can be maintained after the first reionization, thanks to the long hydrogen recombination 
time. Let us make a simple estimate by computing the ratio of the recombination time to 
the dynamical time of a region of overdensity S with ionization fraction x. We find 



for T — 5000 K. It is thus clear that an ionization fraction of order 10~^ — 10~^ will be 
maintained for 5 as high as 10^ (the core density is about 5 — 10^ — 10^), as long as the gas 
is ionized at least to that level at a higher redshift (of course, the first reionization raises 
the gas ionization to a much higher level than this). Such a level of ionization is adequate 

to enable the formation of a sufficient amount of H2 molecules to supply ample cooling at 
T < 10"^K (HAR; Tegmark et al. 1997) and star formation should continue to take place 
in any gaseous halo. Looking forward to the results, it turns out that the condition for H2 
cooling has a rather small effect for the second cosmological reionization, since the amount 
of intergalactic gas able to accrete onto minihalos is very small. 

The halo mass within the virial radius My is related to the virial temperature at z by 
the standard formula: 



where Ty is the virial temperature in Kelvin and Qm is the matter density at zero redshift 
in units of the closure density. Figure 5 relates the virial temperature at each redshift to 
the halo mass for the five indicated virial temperatures. Figure 6 shows the fraction of gas 
in halos as a function of redshift, with virial temperatures greater than the five indicated 
values, = (200, 1000, 2000, 4000, 8000)K. 

The clumping factor in each region affects local atomic processes and is defined to be 
Cj = V'jC'haio + (1 ~ i^i), where Chaio is the effective clumping factor of virialized gas and 
the remaining gas is assumed to have a clumping factor equal to unity for simplicity. We 

will adjust Chaio such that the overall < C > averaged over all regions matches detailed 
simulations at some appropriate redshift. Note that the adjustment on Chaio is necessary in 
order to take into account complicated radiation shielding effect affecting halos. 

Having defined the relevant formulae let us go back to examine the main evolution equa- 
tions more closely. Equation (12) is the rate equation that mass in region i is photoionized; 
i.e., the rate of the amount of gas swept by the ionization front. The sum of ionized gas over 
all the phase space regions at each time step creates a separate (HII) region in the phase 




(19) 




(20) 



- 21 - 




Fig. 5. — relates the virial temperature to the mass of the halo as a function of redshift, for 
five indicated virial temperatures, % = (200, 1000, 2000, 4000, 8000)K. 
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Fig. 6. — shows the fraction of gas in halos as a function of redshift, for five indicated virial 
temperatures, T„ = (200, 1000, 2000, 4000, 8000)K. 
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space with a mean temperature equal to 

N 

Tn+1 — 

1=1 

where the first term inside the parentheses on the right hand side is the temperature of 
region i and the second term accounts for photo-heating, with T^on set to 2 x 10*^K (varying 
Tion from 1 X 10"^K to 3 x lO^'K does not significantly change the results) We set the initial 
neutral fraction of the newly created Hll region to Vn+i = 10^^ , although results obtained 
do not depend on the exact choice for this value, as long as it is a small number. The initial 
mass fraction of this new HII region is equal to 

fN+i = ^^_At, (22) 
at ny 

where At is the time step. We formulate that all phase space regions co-exist in real space 
when averaged over a sufficiently large volume and photoionization ionizes all phase space 
regions proportionally, as indicated by Equation (12). Figures (13,14) are the energy equation 
and rate equation for neutral hydrogen in each phase space region. 

The attentive reader may have noticed the central idea behind this simple approach. 
That is, this approach attempts to mimick the stochastic percolation process during cosmo- 
logical reionization, where individual Hll regions are created and their subsequent evolution 
followed. New HII regions created may originate from a combination of regions, some of 
which have previously been ionized but subsequently recombined, some of which were never 
ionized and some of which may be recently created HII regions but now get percolated. 
We do not include a photoionization term on the right hand side of Equation (14), rather, 
equivalently, we remove the fraction of gas in that region that is being ionized and add it to 
a new phase space region according to Equations (12 and 22). This is equivalent to saying 
that gas is ionized when ionization front passes through, and is not ionized uniformly. 

One important difference between this method and another commonly used method 
which computes the filling factor of HII regions (e.g., Shapiro & Giroux 1987; Haiman & 
Loeb 1997; Madau et al. 1999; Miralda-Escude, Haehnelt, & Rees 2000) is that we follow 
the evolution of all regions including HI, HII regions and partially ionized regions. This 
is important for our purpose, chiefiy because at high redshift HII regions experience rapid 
recombination and cooling. Needless to say, this method is still highly simplified. More 
sophisticated calculations would require simulations that follow detailed three-dimensional 
radiative transfer (e.g., Norman, Paschos, & Abel 1998; Abel, Norman, & Madau 1999; 
Razoumov & Scott 1999; Kessel-Deynet & Burkert 2000; Gnedin & Abel 2001; Gen 2002) 
but will be carried out in our future work. 
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We start a simulation at a very early redshift, say, z — 100, when the gas may be safely 
assumed to cold and neutral. For simplicity we set the initial gas temperature to be equal 
to the microwave temperature, although the gas may be cooler than that, but the results 
would not change noticeably. Thus, the initial condition is this: Ti = 2.73(1 + -2j)K (where 
Zi is the starting redshift), t/i = 1, Ci = 1, /i = 1 and the initial number of entries in the 
phase space (T, y) is iV = 1. N increases linearly with the number of time steps taken. We 
can afford to make the timesteps sufficiently small to be able to follow the ionization and 
energy equations accurately in this simple approach. 

To sclf-consistcntly compute the IMF, wc follow the metal enrichment history of the 
intergalactic medium. Most importantly, there is a critical transition at some point in 
redshift in the metallicity of the IGM, when a certain amount of Pop III stars have formed. 
Subsequent to that transition, sufficient metal cooling would cool collapsing gas inside halos 
to temperatures lower than achievable by H2 cooling, resulting in the "normal" Pop II stars 
with a "normal" Salpeter-like IMF. According to Oh et al. (2001) this transition occurs, 
when a fraction, 3 x 10~^ — 1.2 x 10^^, of the total gas is formed into VMS with masses 
in the range 140 — 260 M0. Without fine tuning, it is likely that VMS with masses outside 
this range should exist. Bearing that in mind we assume conservatively for our calculations 
a threshold fraction of 1 x lO"'*. We note that between the formation of stars and metal 
enrichment of the IGM there should be a delay in time. To accommodate this time lag we 
assume ^fltn (where tH is the Hubble time at that time) after a fraction 1 x 10~^ of total gas 
has formed into VMS, the transition from Pop III VMS to normal stars with Salpeter-hke 
IMF takes place. 

The emission spectrum of metal-free VMS with M > 100 M© is relatively simple and 
has, to a very good approximation, a blackbody spectrum with an effective temperature 
of ~ 10^-^ K, radiating at the Eddington luminosity LEdd = 1-3 x 10^^(M/ Mo)erg s~^ 
for about 3 x 10^ yr, which translates to a hydrogen ionizing photon production rate of 
1.6 X lO^s photons s"^ Mq"^ for VMS (El Eid, Frickc, & Ober 1983; Bond, Arnett, & 
Carr 1984; Bromm, Kudritzki, & Locb 2001; Schaerer 2002). The hydrogen ionizing photon 
rate for normal Pop 11 stars with Salpcter-likc IMF with a low metallicity is assumed to 
be 8.9 X 10"^^ photons s~^ Mq"^, corresponding to a UV emission efficiency of 1 x 10""^ (in 
units of the rest mass energy of the total amount of stars formed) , adopted from Bruzual- 
Charlot stellar synthesis hbrary (Bruzual 2000) for metal poor stars. With the high effective 
temperature VMS are efficient emitters of ionizing photons with approximately ten times 
more hydrogen ionizing photons per unit stellar mass than stars with the Salpeter IMF. 

We use a very conservative star formation efficiency for minihalos with H2 cooling 
c^2 = 0.002 (Abel et al. 1997, 2000; Bromm et al. 2001,2002), the fraction of gas formed 
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into stars out of an amount of gas virialized in minihalos. Star formation efficiency in fialos 
with efficient atomic cooling is assumed to be cj^j = 0.1. 

4.2. Normalizing Ionizing Photon Escape Fraction 

The most important factors determining the ionization process are c*, /eg, C. The first 
two can be combined (see Equation 15). Thus, effectively, there are two primary factors, 
c*/es and C. For the sake of explanation and taking a very conservative stance we use a ratio 
Chi/ch2 — 50, as indicated in the previous subsection, to maximize the radiation emission for 
the second rcionization when Pop II stars are responsible and minimize the emission from 
Pop III stars that arc responsible for the first rcionization. Consequently, we are left with 
the freedom to adjust /es and C. 

Both theoretically and observationally, we have essentially no direct knowledge about 
/es at high redshift in question. One may argue that smaller galaxies at high redshift would 
result in higher escape fraction. On the other hand, one may also argue that higher density 
at high redshift may yield lower escape fraction (e.g., Ricotti & Shull 2000). It thus seems 
most productive to seek empirical constraints on /eg. Our approach is to require that the 
universe is reionized at 2; = 6, as observations suggest (e.g.. Fan et al. 2001; Becker et al. 
2001, Barkana 2001; Cen & McDonald 2002; Litz et al. 2002). With such a normahzation 
point and the fact that the clumping factor C is well constrained by Lya forest observations, 
we are able to tightly constrain the range of /eg. Since the same range of waves in the density 
fiuctuation power spectrum traced by the Lya forest are mostly responsible for the clumping 
of the IGM at ^ = 6. This direct observational constraint is very powerful. 

Figure 7 shows the required ionizing photon escape fraction from galaxies as a function 
of the clumping factor of the IGM at 2; = 6. The cross-shaded area indicates the range 
of the clumping factor constrained by Lya forest observations (Croft et al. 2002), which is 
obtained as follows. Based on the latest high signal-to-noise Lya forest observations. Croft 
et al. (2002) give a constraint on cig = 0.78 ± 0.22 (2cr) for a flat cold dark matter model 
with VLm = 0.25 and A = 1 — Qm- Using slightly different values of Qm does not significantly 
change the range of erg. We use simulations to empirically obtain the clumping factor of HIT 
regions at 2; = 6. The high resolution simulations of Gnedin & Ostriker (1997) are the best 
simulations for this purpose. They give a clumping factor for HIT regions Chu = 35 at 2; = 6 
for a low (J8 ACDM model with os = 0.67, Qm = 0.35, h = 0.70 (see Figure 2 of Gnedin 
& Ostriker 1997). Since the ACDM model at the high redshift behaves like an Einstein-de 
Sitter universe and as a result we can scale their results approximately as follows. We can 
obtain the clumping factor at 2; = 6 for a erg = 0.56 model by using the clumping factor 
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Fig. 7. — shows the required ionizing photon escape fraction from galaxies as a function of 
the clumping factor of the IGM at z = 6. The cross-shaded area indicates the range of the 
clumping factor constrained by Lja forest observations (2cr; Croft et al. 2002). 
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Bit z — 7.4 of the Gnedin & Ostriker (1997) simulation and similarly the clumping factor 
at 2; = 6 for a (jg = 1.0 model by using the clumping factor at z = 4.0 of their simulation. 
Constraining to have an ionized universe at 2; = 6 wc obtain Chii = 24 — 60 at z = 6 for 
the ACDM model with (Tg = 0.56 — 1.0 {2a), shown as the cross-shaded region in Figure 7. 
The solid dots in Figure 7 are obtained by running simulations of the reionization process 
using the method described §4.1 by requiring that reionization complete at 2; = 6 and the 
clumping factor at 2; = 6 matches the value on the x-axis. 

Figure 8 shows the epoch of the first reionization as a function of the normalized escape 
fraction. The cross-shaded region is constrained by the normalization that the universe is 
reionized at 2; = 6 (see Figure 7). We see that, the required constraint by the Lya forest 
observations limits the first reionization epoch to 2; = 15 — 16. 

We note that the required ionizing photon escape fraction of 0.15 — 0.23 from galaxies 
at 2; ~ 6 is significantly higher than those of local starburst galaxies. It should stressed that 
the clumping factor predicted by Gnedin- Ostriker simulations, though high resolution, is 
nevertheless still underestimated. Thus, the actual required ionizing photon escape fraction 
from galaxies may be somewhat higher than indicated. For local starburst galaxies Hurwitz 
et al. (1997) give U < 0-032, 0.052, 0.11 {2a) for Mrk 496, Mrk 1267, and IRAS 08339+6517 
(< 0.57 in the case of Mrk 66). Dcharvcng et al. (2001) give an escape fraction of fcs < 0.062 
for Mrk 54. Heckman et al. (2001) give f^s < 0.06 from local starburst galaxies. Theoretical 
models (Dove & Shull 1994; Dove et al. 2000) give an estimate of fcs = 0.02 — 0.10. Thus, we 
conclude that, in the context of ACDM model, galaxies at 2; ~ 6 appear to demand a higher 
ionizing photon escape fraction than local starburst galaxies. Upon extrapolation it seems 
plausible that a higher ionizing photon escape fraction may be expected for galaxies at still 
higher redshift, those responsible for the first cosmological reionization. For simplicity and 
being conservative we have assumed that the ionizing photon escape fraction from galaxies 
at 2; > 6 is equal to that as determined in Figure 7. 



4.3. Detailed Evolution of Intergalactic Medium 

Let us now examine the evolution of IGM in the context of standard cold dark matter 
model. We use the following model parameters: Cjjj = 0.002, Cjj = 0.1, < C >= 25 at 2; = 6, 
fes — 0.15. This model lies at the lower (left) bound of the cross-shaded region in Figure 7. 
Any model to the right of this model would require higher ionizing photon fraction and thus 
would yield a higher first reionization redshift. A spatially fiat cold dark matter cosmological 
model with = 0.25, fi^ = 0.04, A = 0.75, Ho = 72km/s/Mpc and as = 0.8 is used. The 
results do not sensitively depend on small variations on the cosmological model parameters. 
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Fig. 8. — shows the epoch of the first reionization as a function of the normahzed ioniz- 
ing photon escape fraction from galaxies. The cross-shaded region is constrained by the 
normahzation that the universe is reionized at z = 6 (see Figure 7). 
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We remind the reader that in our formahsm described in §4.1 we have the freedom to adjust 
C'haio to match the required < C > at 2; = 6; we find that Chaio = 702 with the adopted 
cosmological model provides a match to the required < C >— 25 at 2; = 6. 

Figure 9 shows the global hydrogen neutral fraction and the complimentary ionized 
fraction as a function of redshift. The first reionization at z = 15.5 as well as the sustained 
ionized state until z = 13.2 is made possible by Pop III VMS. The redshift z = 13.2 marks 
the transition from Pop III stars to Pop II stars, causing the emission of hydrogen ionizing 
photons to plunge. The suddenly reduced hydrogen ionizing photon emission is no longer able 
to counter the rapid hydrogen recombination process, resulting in the second cosmological 
recombination aX z = 13.2. Since a very small amount of neutral fraction suffices to blank 
out all LycK emission, the universe essentially becomes opaque to Lya photons from z — 13.2 
to 2; = 6. 

However, hydrogen is significantly ionized with nmi/'n-H > 0.6 throughout the long 
second reionization process from z — 13.2 to 2; = 6. To facilitate a better understanding, it 
is useful to show some important time scales involved. We show in Figure 10 the ratio of 
hydrogen recombination time over the Hubble time and the ratio of Compton cooling time 
over the Hubble time. It is noted that at the redshift and density of interest Compton cooling 
dominates over other cooling terms, although adiabatic cooling starts to become important 
approaching the end of the second reionization period. We immediately see that the hydrogen 
recombination time is significantly longer than the Compton cooling time, both of which are 
shorter than the Hubble time at 2; > 8. Therefore, the ICM at early times heated up by the 
photoionization would subsequently cool down more rapidly than rccombining, resulting in 
overcooled but significantly ionized ICM, in the absence of subsequent photoheating. 

More results are shown in Figures (11,12,13,14) for the evolution of mean temperature, 
clumping factor, stellar fraction and number of ionizing photons, respectively, as a function 
of redshift. In Figure (11) we see a sustained ascent in mean temperature of the ICM 
from 2; ~ 30 up to the redshift of the first reionization z — 15.5. Subsequently, the mean 
temperature of the ICM is roughly maintained at lO^K, due to the works of two counter- 
balancing effects: substantial coohng reduces the Jeans mass of the IGM and would increase 
the star formation rate, which would then provide increased heating by photoionization. 
Thus, the mean temperature of the IGM is self-regulated due to the competition between 
the Compton cooling and photoheating. For the cold dark matter cosmological model it 
happens that the star formation activities are at a level such that the two balancing terms are 
comparable in magnitude, resulting in a fairly mild evolution of the IGM mean temperature 
seen in Figure 11. 

From Figure 12 it is evident that in the period preceding the first reionization from 
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Fig. 9. — shows the global mean of the hydrogen neutral (solid) and complimentary ionized 
(dashed) fraction as a function of redshift. 
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Fig. 10. — shows the ratio of the recombination time to the Hubble time (sohd curve) and 
ratio of the Compton coohng time to the Hubble time (dashed curve) a function of redshift. 



- 32 - 




Fig. 11. — shows the evolution of the mean IGM temperature as a function of redshift 
during the second cosmological reionization process. The sohd vertical tick indicates the 
first reionization epoch. The dashed vertical tick indicates the transition epoch from Pop III 
stars to Pop II stars. The dotted vertical tick indicates the second reionization epoch. 
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Fig. 12. — shows the evolution of the mean clumping factor as a function of redshift during 
the second cosmological reionization process, consistent with direct numerical simulation of 
Gnedin & Ostriker (1997). The solid vertical tick indicates the first reionization epoch. The 
dashed vertical tick indicates the transition epoch from Pop 111 stars to Pop 11 stars. The 
dotted vertical tick indicates the second reionization epoch. 
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Fig. 13. — The left y-axis shows the evolution of the fraction of gas formed into stars as a 
function of redshift. The right y-axis translates the fraction of gas formed into stars into 
the number of ionizing photons per baryon as a function of redshift. The solid vertical tick 
indicates the first reionization epoch. The dashed vertical tick indicates the transition epoch 
from Pop III stars to Pop II stars. The dotted vertical tick indicates the second reionization 
epoch. 
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Fig. 14. — shows the cumulative number of ionizing photons per hydrogen atom produced 
as a function of redshift. The sohd vertical tick indicates the first reionization epoch. The 
dashed vertical tick indicates the transition epoch from Pop III stars to Pop II stars. The 
dotted vertical tick indicates the second reionization epoch. 
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z — 18 to z — 15.5 a large fraction of gas is unable to accrete onto minihalos due to heated 
temperature of the IGM and results in a large decrease in the clumping factor of the IGM, 
followed by a mild increase in the clumping factor in the brief ionized state from z = 15.5 to 
z = 13.2. From z = 13.2 to z = 6 the clumping factor of the IGM increases steadily, parallel 
to the increase in the nonlinear mass scale with time, which is then followed a very brief, 
moderate drop immediately after z = 6 when there is an increase in the temperature of the 
IGM from T ~ lO^'K to T ~ 2 x lO^'K (see Figure 11) due to the completion of the second 
reionization. 

Some interesting features may be found in Figures (13,14). We see in Figure 13 that, 
although stars formation in minihalos (dashed curve) dominate over that in large halos (solid 
curve) at, z > 17, the latter gradually takes over and become dominant at, z < 17. Therefore, 
while Pop III stars in minihalos may be a very important source of heavy elements to enrich 
the IGM, the first cosmological reionization is largely due to Pop III stars in large halos. 
As a result, at least for this particular case, the first reionization would have occurred even 
without the validation presented in §2,3 that H2 cooling is sufficient in minihalos during the 
first reionization process. About 6 x 10^^ of the total gas has formed into Pop 111 VMS by 
the time the universe was first reionized. We note that we have adopted a very conservative 
ratio of c^i/c^2 ~ 0.1/0.002 = 50 for this example. If, for example, the ratio were 10, then 
Pop III stars in minihalos would largely be responsible for first reionization. By the time the 
universe was reionized for the second time, about 1% of the total gas has formed into stars. 
From Figure 14 we see that at the epoch of the first reionization about 4 ionizing photons 
per baryon have been produced by Pop III stars. About 10~^ fraction of gas formed into 
Pop III stars at a slightly later epoch and the delayed transition from Pop III stars to Pop II 
stars occur at z = 13.2, clearly seen in Figures (9,11,12). More than 60 ionizing photons per 
baryon have been produced by the time than the universe is reionized for the second time 
at z = 6. 

In summary, the cosmological reionization process has several interesting characteristics. 
First, the universe was generically reionized twice, one at 2; = 15 — 16 by Pop 111 VMS and 
the other at 2; ~ 6 by normal Pop II stars. Most of the Pop 111 VMS responsible for the 
first cosmological reionization would come from large halos if star formation efficiency in 
minihalos is less than a few percent of that of large halos; in this case, our argument that H2 
cooling is efficient in minihalos during the first cosmological reionization, while still true, may 
not be necessary. Otherwise, most of the Pop III VMS responsible for the first cosmological 
reionization would come from minihalos, where efficient H2 cooling is required. Second, the 
IGM is maintained at a rather "warm" temperature T ~ 10^ K from the first reionization 
through the second reionization; i.e., the universe is significantly heated up even in this "dark 
age" preceding the second, complete reionization of the universe at ^ ~ 6. Third, the IGM 
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is kept at a significantly ionized state throughout the second reionization period. Finally, 
the overwhelming fraction of star formation activity responsible for producing the ionizing 
photons in the period from the first reionization through the second reionization are in large 
halos. 



5. Discussion 

The possible indication of a more top-heavy initial stellar mass function (IMF) in early 

galaxies or earlier stages of galaxies than prcscnt-day IMF (Salpeter 1955) was suggested 
five decades ago (Schwarzschild & Spitzer 1953). Cosmological consequences of the first 
generation, massive stars at high redshift have been discussed in a variety of contexts (e.g., 
Layzer & Hively 1973; Carr 1977; Rccs 1978; Rowan-Robinson, Ncgroponte, & Silk 1979; 
Puget & Heyvaerts 1980; Tarbet & Rowan- Robinson 1982; Carr, Bond, & Arnett 1984; 
Haiman k Loeb 1997; Barkana & Loeb 2001). 

The reionization picture presented here would have a wide range of profound implica- 
tions for many aspects of structure formation. There are many questions that need to be 
addressed. It is beyond the scope of this paper to explore any of these issues in significant 
details and we will study these and other relevant issues in subsequent investigations. But 
we will give some simple estimates or analyses for some selected issues. 

5.1. Initial Metal Enrichment of the Intergalactic Medium 

Direct observational evidence for massive Pop III stars may be just emerging recently. 
While Luck & Bond (1985) and others have previously indicated the need for VMS (M > 
100 Mq) to explain the overabundant a elements, Qian, Wasserburg and collaborators (Wasser- 
burg & Qian 2000; Qian & Wasserburg 2000,2002; Oh et al. 2001) have recently stressed the 
unique signature of Pop III VMS and suggested that Pop III stars could promptly produce 
the observed abundance patterns of metal poor stars (McWilliam et al. 1995; Ryan, Norris, 
& Beers 1996; Rossi, Beers, & Sneden 1999; Norris, Ryan, & Beers 1997,2001; Norris et al. 
2002; Hill et al. 2002; Depagne et al. 2002). 

Given the shallow potential wells of Pop III galaxies, supernova explosions would likely 

blow away the ejectas along with a large fraction of the diffuse interstellar gas (Mac Low 
& Ferrara 1999; Mori, Ferrara & Madau 2002). From Figures (9,11-14) we see that the 
transition from Pop III stars to Pop II stars occurs a later time than the first reionization; 
i.e.. Pop III stars ionized the universe for the first time. Oh et al. (2001) have suggested 
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this possibility based on the observed transition of abundance pattern at [Fe/H]=-3 in metal 

poor stars summarized by Qian & Wasserburg (Wasserburg & Qian 2000; Qian & Wasserburg 
2000,2002), who found that at lower value VMS dominates the enrichment, whereas at higher 
value a sharp rise in the abundances of the heavy r-process elements such as Ba and Eu in 
galactic halo stars with [Fe/H]> —3 signifies the occurrence of Type II supernovae of normal 
stars with masses of 10 — 60 Mq. 

One may expect to see the Pop III star abundance patterns in low density regions of 
the universe, such as Lja clouds or voids, where subsequent additional enrichment may be 
minute. Recent work on Lja clouds suggests that metal enrichment by Pop III stars be 
consistent with Lya cloud observations and the necessary enrichment occur prior to ^ ~ 4.6 
(the highest epoch for the available Lya cloud data analyzed) (Qian, Sargent, & Wasserburg 
2002). 

One obvious advantage for Pop III supernovae to enrich the IGM with the first metals 
at ;2 3> 4 is that it is much easier to relatively uniformly disperse the metals across the IGM 
for two reasons. First, the disturbances to the density and velocity fields are smaller at high 
redshift, because each Pop III galaxy needs to fill only a very small IGM volume. Second, 
any significant large-scale motions would decay away rapidly, in the absence of dynamical 
support. As a result, the excellent agreement found between simulations and observations 
with respect to the Lja forest at 2; = 2 — 4 (Ccn ct al. 1994; Zhang et al. 1995; Hernquist ct 
al. 1996; Miralda-Escudc ct al. 1996; Bond & Wadsley 1997; Thcuns et al. 1998) would not 
be altered, although recent simulations also indicate that low redshift galactic winds would 
not spoil the main properties of Lya forest produced by previous simulations (Theuns et al. 
2002). 



5.2. Intergalactic Magnetic Field 

Another possible consequence of supernova explosion and expulsion of the gas into the 
IGM is the magnetic feedback to the IGM. Rccs (1994) has first pointed out the impor- 
tance of resident supernova remnants in galaxies being a substantial large-scale seed field 
for galactic dynamos. Let us give a simple estimate here for the contribution of the Pop 
III supernova remnants to the intergalactic magnetic field. FoUoing Rees (1994) we should 
use local observations as a guide. The Crab Nebula (a plerion) has a magnetic field of 
strength B ~ lOO/iG currently occupying a sphere of radius r ~ 1 pc, totaling a flux of 
V'Crab ~ 3 X 10^^ G cm^. Assuming flux conservation and Pop III magnetic bubble filling 
factor of unity give 5igm = (^^^^nS?^)'^'( Ji&)^Crab, where Mui ~ 100 Mq is the mass 
of a typical Pop III star, Mcrab ~ 5 Mq is the zero-age main sequence progenitor star for 
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Crab; /m ~ 2.5 x 10~^ is the fraction of baryons formed into Pop III stars (see Figure 14) by 
redshift z = 13.2 when the transition from Pop III to Pop II occurs; pi, is the mean baryonic 
density at zero redshift; Zgj is the redshift of ejection of the magnetic field into the IGM; 
we have assumed that the magnetic flux is approximately proportional to the stellar mass. 
Inserting all the numbers gives Biqm{z = 13.2) ~ 1 x 10~^ G at z = 13.2 and subsequently 
Bigm{z) ~ 1 X 10~^(U|)^ G. This magnetic field, having a significantly larger amplitude 
than that produced by gravitational shocks in the collapse of large-scale structure (Kulsrud 
et al. 1997), could serve as a seed field for subsequent galaxy formation. The mean separa- 
tion between Pop III galaxies is of order 100 comoving kpc. Thus, if Pop III supernovae are 
responsible for enriching the IGM relatively uniformly to a metallicity of about a thousandth 
of the solar value, it is also likely that the magnetic field lines originating from supernovae 
remnants would be stretched to fill up the intcrgalactic space, resulting in an initial magnetic 
field possibly coherent on scales as large as ~ lOOkpc. 

In addition, miniquasars powered by Pop III black holes may produce mini radio jets, as 
would have been implied by the observational fact that radio jets are observed in accretion 
disks on a wide range of scales in a wide variety of astrophysical systems. The magnetic 
field from miniquasars may be as important as or more important than that from the Pop 
III stars; but a reliable estimate is difficult. 

5.3. Pop III Black Holes 

Without fine tuning the IMF of Pop III stars, it seems likely that ~ 200 Mq black hole 
from Pop III stars more massive than 260 Mq and smaller black holes (Mbh ~ 10 — 50 Mq) 
from Pop III stars less massive than 140 Mq would form. The possible bimodality of the 
distribution of black hole mass is interesting but the consequences are too complicated to be 
easily outlined. Those Pop III black holes will be building blocks for subsequent structures 
and the question is: where will they go? In general, since there is no correlation between 
the small-scale structures that form Pop III galaxies and later structures, one would expect 
to find those black holes in all environments, including globular clusters, galactic disks, 
galaxy halos and intergalactic space. It is possible that the dynamic formation of halos and 
galaxies in hierarchical structure formation may be significantly altered in the presence of 
these massive black holes, especially in the cores of the relevant structures, such as globular 
clusters, galaxies, etc. The number density of Pop III black holes of mass ~ 100 Mq is of 
order 10^ - 10^ per comoving Mpc^, comparable to the number density of globular clusters 
and mass density comparable to present mass density in supermassive black holes in the 
centers of galaxies (e.g., Merritt & Ferrarese 2001). The interesting agreement between this 
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result and that of Madau & Rees (2001) with regard to the mass density of Pop III black 
holes is traceable to the fact that the 3a density peaks used in Madau & Rees's (2001) model 
happen to yield about the same collapsed fraction as ours. From the inferred number density 
of Pop 111 black holes it is quite possible that many globular clusters could be seeded by 
those black holes. There is clearly no shortage of Pop 111 black holes to provide seeds for 
later supermassive black holes seen in centers of local galaxies (e.g., Tremaine et al. 2002), 
AGN and quasars (e.g., Rees 1984,1990). 

How could these black holes be observed? Miralda-Escude & Gould (2000) pointed out 
the possible existence of stellar black hole clusters at the Galactic center. Their reasoning 
may be equally applied here. Since Pop 111 star mass fraction is of order of 10~^ of total 
baryonic mass, as our calculation indicates, their should be ~ 10^ black holes either having 
mass 100 or 10 M© (noting that the mass fraction in black holes will be a factor about 
10 lower for less massive black holes). The reader is encouraged to refer to Miralda-Escude 
& Gould (2000) and also Sigurdsson & Rees (1997) for a detailed discussion of possible 
manifestations of such a black cluster, bearing in mind that the black holes here are more 
massive and may enhance some of the effects discussed. 

In the hierarchical structure formation process, these Pop III black holes are expected 
to grow (there are a large number of e-folding times available) and merge, along with larger 
cosmic structures. The merger of black holes may be detectable by the gravitational wave 
experiment Laser Interferometer Space Antenna (LISA) (e.g., Menou, Haiman, & Narayanan 
2001; Hughes 2002). 



5.4. Detecting Pop III Galaxies and Pop III Hypernovae 

We will examine the observability of Pop 111 galaxies and their associated hypernovae 
with respect to SIRTF and JWST (James Webb Space Telescope). 

The physical sizes of the galaxies in the redshift range 2; = 13 — 15 is unknown but it 
seems unlikely that they will be larger than their lower redshift counterparts. Steidel et al. 

(1996) find that the great majority of 2; > 3 objects have a half-light radius ~ 2 kpc, which, 
if placed at Zs = 17, would subtend ~ ".05 (for the adopted model), more than a factor 
of ~ 10 below the angular resolution of SIRTF. In other words, these high redshift galaxies 
are point sources to SIRTF. The flux density (the power per unit antenna area and per unit 
frequency interval) of a point source (without cosmological surface brightness dimming) at 
frequency v is (Weinberg 1972): S = where zs is the source redshift and u is the 

observed photon frequency. P is the intrinsic power, the power emitted per unit solid angle 
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and per unit frequency interval, related to the luminosity L of the source by P = L/AttAu, 
where Ai/ is the rest frame bandwidth at which the source has a luminosity L; dcm is the 
comoving distance to the source. Using a rest frame band 0.09 — 0.27//m (corresponding to 
the observer's band 1.3 — 3.8 iim) and = 13-2 we obtain 

S = 1-45(y^) nJy. (23) 

The point source detection sensitivity at la level (for relatively long integration time > 
500 seconds) of the Infrared Array Camera (IRAC) on SIRTF at 3.5 — 4.5 fim is (SIRTF 
Project 1997) 6*3.5 ~ 1/ (100/t) + 0.6^ /xJy, where t is the integration time in seconds. The 
second term inside the square root is the confusion noise limit due to faint unresolved sources, 
modeled by Pranceschini et al. (1991). 

The total luminosity of a Pop 111 galaxy ranges from L = Ix 10'^(42/0.002) (Mh/10^ Mq) L, 
for minihalos to L = 1 x 10^°(cHj/0.1)(Mh/2 x 10^ M©) L© for large halos; only a small frac- 
tion of the total luminosity is in the indicated band. It is clear that no Pop III galaxies hosted 
by minihalos will be directly detectable by SIRFT, even without the confusion noise hmit. 
Perhaps only the high mass end of Pop III galaxies with mass > 10^ Mq (which would give 
> 0.7//Jy) may be detectable by SIRTF, in the presence of confusion sources. The sensitivity 
of the JWST is close to InJy (Stockman & Mather 2000). Therefore, Pop III galaxies with 
large halos will be detectable by JWST, whereas Pop 111 galaxies with minihalos may still 
be beyond reach of even JWST. Our calculations indicate a large fraction of Pop 111 galaxies 
will be detectable by JWST; a more detailed treatment will be reserved for a future paper. 

The total number of Pop 111 hypernovae (HN) is Nun = 2.5 x 10-^pb{An/3)Rfi/Mui = 
1.3 X 10^^ where Mm = 100 Mq for Pop 111 star mass is used, Rh = 6000/^"^ Mpc is the 
comoving Hubble radius and pf, is the comoving mean baryonic density. We can estimate 
the Pop 111 hypernova surface number density at any given time to be 

Ehn = 1.3 X 10^6 X ( ( L_ ] X 0.60 ^ 2.3 arcmin-^ (24) 

V y V47r X 1.2 X lOWminV ' 

where A^hn is the time duration of each hypernova event in the hypernova restframe and 
in is the age of universe at z — 13.2; the first parenthesized term in the above equation 
takes account the finite duration of each supernova event and an intrinsic duration for the 
hypernovae A^hn = 1 yr is used; the last term 0.60 takes into account that the observable 
redshift interval z > 13.2 is 60% of the total volume of the universe then. This HN surface 
density is comparable to the surface density of SNe derived by Miralda-Escude & Rees (1997), 
based on the observed metallicity in the Lya clouds. 

Assuming the optical luminosity of a hypernova to be 1 x 10^° L© (about 10 times that 
of the local normal supernovae), the fiux density of Pop HI hypernovae at 2; = 13.2 will 
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be at a level ~ 0.014//Jy. Even without background confusion {if the variable source has a 

variability time scale shorter than the indicated integration time), in order to detect such a 
flux density at la level, an integration time of 6 days on SIRTF IRAC will be required, which 
seems impractical. Even if this long integration is achievable, the fact that the variability 
time scale of a hypernovae at the considered redshift may be quite long indicates that it 
will be extremely difficult to directly detect high redshift hypernovae, unless their flux is 
significantly amplified. 

In Appendix A we derive the probability of strong gravitational lensing by massive 
clusters of galaxies and find that a random source behind a massive galaxy cluster will have 
a probability 

6.7 X 10-2 / Qpov 



PciuM = ——2 ^^^^^ (25) 

/x^ V25 arcmm / 

being magnified by > where /i is the lensing magnification and Jlpov is the field of view of 
a telescope. The number of Pop III hypernovae, magnified by > /i, per field of view centered 
on a massive cluster is then 

NhN = -Pclust(A*)f^FOvSHN- (26) 

In order for the IRAC camera of SIRTF to detect at na statistical level a point source, 
which would have a fiux density of S in the absence of gravitational lensing magnification, 
one requires that 

Sn = ny/TooJi, (27) 

where point source confusion limit is removed; t is integration time in seconds. We plot 
the number of hypernovae, A^HN,det) per field of view detected at 3(7 confidence level against 
integration time t in Figure 13, assuming Lhn = 1 x 10^° Lq giving unlensed flux density of 
S = 14.5nJy at the observed wavelength of several fim. To put the matter in perspective, 
1,000 target fields each with lOk-second integration time would detect ~ 300 multiply lensed 
Pop III hypernovae at 3a confidence level. 

We note that since only highly magnified Pop 111 hypernovae will be observable, the 
observed candidates should have distinct features: the multiple images will significantly help 
the identification process. If the cluster lens is known, its lens potential derived independently 
elsewhere would further constrain the image configurations. Once such images are found, 
periodic monitoring up to a few years would eventually verify the transient nature of the 
Pop III hypernovae. Since the total number of detectable hypernovae is just proportional to 
the total integration time, one could choose an optimal strategy such that the observation 
is most sensitive to the expected splittings of the images: a shorter exposure will require a 
higher magnification from a more central region of the lens, where splittings may decrease 
due to a smaller velocity dispersion at the center of the cluster. 
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Fig. 15. — shows the average number of strong gravitational lensing magnified Pop hyper- 
novae that will be observed per IRAC field of view (25 arcmin^) as a function of integration 
time. 



-44- 



At a flux density of a level ~ lOnJy, Pop III hypernovae will be detectable by JWST 
(Stockman & Mather 2000) without gravitational lensing magnification. The possible ob- 
servable duration of 4 — 5yrs for the initial bright phase of Pop III hypernovae, corresponding 
to 80 — 90 d intrinsic duration (Woosley & Weaver 1996) will be a signature. 

5.5. EflFects on the Cosmic Microwave Background 

The earlier reionization would significantly increase the electron Thomson scattering 
optical depth. This is so because the IGM ionization fraction remains appreciable during 
the period between the first reionization {z = 15.5) and the second reionization {z = 6.1). 
Figure 16 shows the cumulative Thomson scattering optical depth as a function of redshift. 
We see that the Thomson scattering optical depth increases from Te{z < 6.1) = 0.027 
to Te{z < 20) = 0.097. Evidently, the first reionization boosts the optical depth of the 
second reionization by a factor of ~ 3, making the total Thomson optical depth much more 
significant. Given the narrow allowed range of the first reionization epoch between z — 15—16 
(Figure 8), our prediction of Tg ~ 0.10±0.03 [where AtaUe — 0.03 is approximately estimated 
considering uncertainties on C{z), fesi^), c*{z) and fib] is fairly robust, unless the ionizing 
photon escape fraction from the Pop III galaxies is extremely small compared to what is 
required at 2; ~ 6 to reionize the universe (Figure 7). The implications for the CMB (for an 
excellent recent review see Hu & Dodelson 2002), in particular, for polarization of the CMB 
(e.g., Seljak 1997; Zaldarriaga 1997; Kamionkowski et al. 1997) is significant. 

We will follow the recent analysis by Kaplinghat et al. (2002) on the detectability of 
polarization of the CMB. Kaplinghat et al. (2002) show that the Microwave Anisotropy 
Probe (MAP) satellite will be able to measure Tg to an accuracy of 0.02 — 0.03 (Ic). Thus 
it appears that MAP will be able to distinguish the two reionization scenarios at about 
2 — 3(7 confidence level: the standard reionization scenario where the universe is fully ionized 
at 2; < 6 and fully neutral at 2; > 6 (giving Tg = 0.027) versus the reionization scenario 
presented here where the universe is fully ionized at 2; < 6 and partially ionized at 2; = 6 — 20 
(giving Tg = 0.097). According to Kaphnghat et al. (2002) Planck surveyor will probably 
be able to probe the detailed reionization history in addition to discriminating between the 
standard reionization model and the current reionization model at a very high confidence 
level. 
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Fig. 16. — shows the cumulative Thomson scattering optical depth as a function of redshift. 
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5.6. Non-detection of Pop III Stars Locally 

There is a long history of search for Pop III stars (Bond 1970,1981; Hills 1982; Bessell 
& Norris 1984,1987; Carr 1987; Cayrel 1996; McWilliam et al. 1995; Ryan, Norris, & Beers 
1996; Rossi, Beers, & Sneden 1999; Beers 2000; Norris, Ryan, & Beers 1997,2001; Norris et 
al. 2002; Hill et al. 2002; Depagne et al. 2002) or primeval galaxies (Partridge 1974; Davis 
& Wilkinson 1974). If Pop HI stars are as massive as suggested, perhaps it is not surprising 
that no single Pop HI star should have been found, while over 100 metal-poor stars with 
—4 < [Fe/H] < —3 have been found (Cayrel 1996), since all Pop HI stars would end either 
as supernovae or black holes. The proposed model, unfortunately, would mark the end of 
search for Pop HI stars in the local universe. 

5.7. Detectability of GalcLxies Beyond z = 6 

It is interesting to note, from Figure 9, that there is a substantial redshift range beyond 
the end redshift of the second reionization z = 6, where the IGM is already substantially 
ionized. For example, fiui/nu ~ 0.1 — 0.2 at z ~ 6.0 — 7.5. This feature of an extended 
redshift interval of high ionization appears to be generic. This result would have interesting 
consequences on the observability of galaxies beyond the second reionization epoch. Taking 
the result at face value, for a galaxy at z = 6.5, the Stromgren sphere produced by the galaxy 
will be a factor of about 2 larger than that with the case where the galaxy is embedded in a 
completely neutral IGM. This would help reconcile the seemingly conflicting observational 
claims: on one hand, the universe appears to be reionized at 2; ~ 6 when ionizing radiation 
background is seen to rise sharply (e.g.. Fan et al. 2001; Becker et al. 2001, Barkana 2001; 
Cen & McDonald 2002); on the other hand, Lya galaxies at redshift as high as z > 6.5 
(Hu et al. 2002; Kodaira et al. 2003) has been detected (see Haiman 2002 for an alternative 
explanation). The optical depth due to the damping wing of neutral hydrogen outside the 
Stromgren is in general reduced by a factor of {nYn/rin)^^'^ (Miralda-Escude & Rees 1998; 
Cen & Haiman 2000; Madau & Rees 2000), which is a factor of 0.08 for nm/nu = 0.15. 
One prediction is that the number of observable Lya galaxies would quickly thin out beyond 
2; ~ 7, due to the combined effect of a rapid increase of neutral fraction beyond 2; ~ 7 (see 
Figure 9) and a decrease of the number of large Lya galaxies. More detailed treatment of 
this important subject will be done in a separate paper. Observational campaign by several 
groups (Hu et al. 2002; Rhoads et al. 2001; Ajiki et al. 2002) to detect high redshift Lya 
galaxies may be able to probe the detailed reionization structure near the end of the second 
cosmological reionization. 



-47- 



5.8. Hydrogen 21cm Line due to Minihalos Prior to the Second Cosmological 

Reionization: A Test 

The possibility of probing the high-z universe with hydrogen 21 cm Une absorption or 
emission has been suggested and investigated in various contexts by many authors (Hogan 
& Rees 1978; Scott & Rees 1990; Subramanian & Padmanabhan 1993; Kumar, Padmanab- 
han, & Subramanian 1995; Bagla, Nath, & Padmanabhan 1997; Madau, Meikien, & Rees 
1997; Shaver et al. 1999; Tozzi et al. 2000; CariUi, Gnedin, & Owen 2002; Ihev et aL 2002; 
Furlanetto & Loeb 2002). The mean, neutral medium expanding with the Hubble flow at 
redshift z would produce an optical depth for the hydrogen 21 cm resonant absorption is 
(Shklovsky 1960) 

(28) 

where H{z) is the Hubble constant at redshift z; n{z) is the mean atomic hydrogen number 
density at z; gi = 1 and (?2 = 3 arc the statistical weights of the lower and upper hyperfine 
levels; u is the frequency of the hydrogen 21cm line; Tgp is the spin temperature of the atomic 
hydrogen; u is the frequency of the 21 cm line; c is the speed of light; k is the Boltzmann's 
constant; A21 is the spontaneous decay rate of the hyperfine transition of atomic hydrogen; 

and Qft are total matter and baryonic matter density at 2; = in units of closure density; 
h — ifo/lOOkm/s/Mpc, where Hq is the Hubble constant at 2; = 0. Clearly, the uniform 
medium would only cause a very modest absorption at 2; ~ 6 — 15. 

Minihalos would cast significantly larger optical depth. For 21cm absorption by atomic 
hydrogen in virialized minihalos we may rewrite Equation (30) approximately as: 

178n(2:)r^ ^2 ^ hu 

Thyz) = ■ rA2i- — -—— (29 j 

9i{92 + 9i) Snu-^kTsp 

where and cr„ are the virial radius and velocity dispersion within the virial radius, respec- 
tively; the factor 178 is the relative overdensity of matter within the virial radius. Using 
relations cr^ = GMy/2r^, = 178(47r/3)p(2;)r^ [where p{z) is the mean total density at z] 
and the definition of Hubble constant, we can transform Equation (31) to 

ni.) = 2VuSfi.) = 0.023(^)-H^)(5M!)-./.(i±£)3/. (30) 

We see that gaseous minihalos with Tgp ~ 10^ — 10^, if they exist, could produce significant 
optical depths at 21 cm at a level of 0.01 — 0.1 (Furlanetto & Loeb 2002; see also CariUi, 
Gnedin & Owen 2002 for a numerical treatment). Note that the spin temperature Tgp would 
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be bracketed by the cosmic microwave background temperature and the kinetic temperature 
(approximately the virial temperature) of a minihalo. 

Purlanetto & Loeb (2002) have pointed out that the 21 cm forest produced by minihalos 
prior to cosmological reionization will be detectable by the next generation of low-frequency 
radio telescopes, such as Low Frequency Array (LOFAR) and Square Kilometer Array (SKA). 
However, this will only be possible, if the minihalos could retain their gas, as they would 
in the conventional reionization scenario. In the reionization scenario presented here, as 
indicated in Figure 14, no significant amount of gas will be able to accumulate in minihalos 
throughout the second reionization period. Consequently, in our scenario, 21 cm forest lines, 
which would otherwise be produced by minihalos without the first reionization, have been 
largely wiped out. It is noted that 21 cm absorption measurements would require bright 
background radio sources to be present at very high redshift. 

Similarly, all hydrogen 21 cm emission line forest due to minihalos, as many as ~ 100 
lines per unit redshift at z ~ 9, as predicted by Iliev et al. (2002) in the standard reionization 
scenario, will be absent in the present reionization scenario. In addition, the 21 cm emission 
spatial and spectral signature at 2; > 6 (Madau et al. 1997; Tozzi et al. 2000) will also be 
significantly altered in the sense that signals due to minihalos will be removed. 

Hydrogen 21 cm observations, either in absorption or emission, will be a definitive test 
to distinguish between the two reionization scenarios. 

5.9. Metal Absorption Lines at 2; > 6: A Probe of the Second Reionization 
Process and Initial Metal Enrichment History 

Oh (2002) pointed out that some metal absorption lines may be observable at high 
redshift. He identified 01 (1302A) and Sill (1260A) fines, longward of the hydrogen Lya fine, 
as the best candidates for tracing neutral hydrogen due to the proximity of their ionizational 
potentials to the hydrogen ionization potential. Since both these metal species trace closely 
the neutral hydrogen density, 01 and Sill absorption forest may provide a way to probe 
the reionization history. Two possible utilizations are noted here. First, both 01 and Sill 
absorption forests are expected to thin out at 2; < 7, which could provide a direct probe of 
the end of the second reionization process. We note that, for IGM which is primarily metal 
enriched by Pop III stars, the overabundances of a-elements such as oxygen and silicon 
relative to the solar value (Heger & Woosley 2002) would greatly increase the optical depth 
of the oxygen 01 and Sill absorption forests, at the apparent low metallicity (defined by 
[Fe/H] ~ 10~^). This would substantially increase the number of observable 01 and Sill 
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absorption lines predicted by Oh (2002). Second, the relative optical depths of 01 and 
Sill absorption lines may shed light on the relative abundance of O to Si, giving useful 
information on synthesis in Pop III stars (Qian & Wassurberg 2000,2002). In addition, 
gradually increasing contribution from Pop II supernovae (type II) to the IGM metals with 
time may also show up in the redshift evolution of the relative optical depths of 01 and Sill 
absorption lines, possibly providing information on the overall metal enrichment history of 
the universe at those early times. 

5.10. Comparisons with Some Other Works 

The author was made aware of an earlier paper by Wyithe & Loeb (2002), who in- 
vestigated the reionization histories due to stars and quasars. They assume that minihalos 
do not form stars at the redshift of interest but systematically explore the 2-d parameter 
space spanned by the ionizing photon escape fraction and the transition redshift from a VMS 
IMF and a normal IMF. They pointed out that, for some restricted region in the parameter 
space, hydrogen and/or helium II may be reionized twice. The primary differences between 
their study and the current study are three-fold. First, we show that H2 coohng thus star 
formation in minihalos takes place unimpeded, as long as gas is able to accrete onto them, 
whereas they do not consider star formation in minihalos. Second, we constrain the ionizing 
photon escape fraction directly using Lja forest observations. As a consequence, we show 
that it seems inevitable that the universe will be reionized first by Pop III stars either in 
minihalos or in large halos and second by Pop II stars. Finally, we use a new computational 
method to follow self-consistently the evolution of the entire IGM, including HII regions, HI 
regions and partially ionized regions. The two studies are, however, complementary in their 
different approaches. 

Venkatesan, Tumhnson, & ShuU (2003) also examined the consequences of Pop III metal- 
free stars on the reionization of the universe. They indicated the possibility of reionizing He 
II twice with an incomplete first attempt of reionizing He II. 

Mackay, Bromm, & Hernquist (2002) presented a three-phase reionization picture. With 
a top-heavy VMS IMF, they pointed out that the first transition from molecular cooling 
phase to atomic cooling phase occurs at 2; ~ 30, when their demand of 10 Lyman- Werner 
photons per baryon is met. The proposed first transition in Mackay et al. (2002) is dictated 
by the lack of H2 cooling in minihalos. We pointed out here positive feedback and as a result 
there is no such transition. Mackay et al. (2002) find that at 2; = 15 — 20 the IGM has been 
enriched to a level of 10~^'^ Zq that the second transition to a normal IMF occurs. Their 
assessment of the transition epoch from Pop III to Pop II is consistent with but at a slightly 
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higher redshift than our calculation. Their universe, however, is not reionized until z ~ 6. It 
is noted in passing that Mackay et al. (2002) pointed out a novel Pop II. 5 stellar population 
associated with cooling shells of Pop III supernovae. In summary, the primary differences 
between this study and that by of Mackay et al. (2002) are two-fold. First, we indicated 
that the positive feedback mechanisms on II2 formation allow star formation to occur in 
minihalos throughout the first reionization period up to 2; ~ 13.2. Second, we constrain the 
ionizing photon escape fraction by the direct Lya forest observations to eliminate this major 
uncertainty. 



5.11. Uncertainties 

Perhaps the most uncertain in the chain of derivation is the IMF of the Pop III stars. 
However, most of the main conclusions, including positive feedback of Pop III star formation 
and two times of cosmological reionization, are likely to hold for any sufficiently top-heavy 
IMF. While the simulations with sub-solar mass resolution (Abel et al. 2002; Bromm et al. 
2002) have clearly shown the formation of Pop III VMS, it is yet unclear what the angu- 
lar momentum transport mechanisms are. It is possible, in principle, that hydrodynamic 
processes, included in the quoted simulations, could transport the unwanted angular mo- 
mentum outward, as the simulators have advocated. If local star formation is a guide, one 
could imagine that very massive Pop III stars may form in binaries or multiples, in which 
case angular momentum removal would be easily achieved. It is also intriguing to note that, 
if we extrapolate the observed ratio of super massive black hole to bulge mass (e.g., Tremaine 
et al. 2002), then one would expect to find a black hole of mass ~ 100 — 1000 Mq at the 
center of each minihalo of mass 10^ — 10^ Mq. Thus, perhaps one should not be too surprised 
that nature has managed to form a compact object (star) of mass ~ 100 Mq formed at the 
center of a minihalo, if the quoted observations are any empirical guide; it may be argued, 
angular momentum-wise, that it is a less stringent task to collapse a higher density gas cloud 
at higher redshift to form a star than a lower density gas cloud at lower redshift to form a 
black hole, as observed. 

We have ignored the possible contribution from quasars to the second reionization pro- 
cess. If the emission rate of ionizing photons from quasars is proportional to the star for- 
mation rate, as indicated by low redshift observations (Boyle & Terlevich 1998; Cavaliere & 
Vittorini 1998), then, to zero-th order, all that is needed is to rcnormalize /es such that the 
second reionization completes at z ^ 6 and results would remain largely unchanged. 
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6. Conclusions 

The conclusions of this paper consist of two sets. 

I) Wc put forth two new mechanisms for generating a high X-ray background during 
the Pop III era. namely, X-ray emission from the cooling energy of Pop III supernova blast 
waves and from miniquasars powered by Pop III black holes. We show, consequently, II2 
formation in the cores of minihalos is significantly induced, more than enough to compensate 
for destruction by Lyman- Werner photons produced by the same Pop III stars. In addition, 
another, perhaps dominant process for producing a large number of H2 molecules in relic HII 
regions created Pop III galaxies, first pointed out RGS, is quantified here. We show that H2 
molecules produced by this process may overwhelm the Lyman- Werner photons produced by 
the stars in the same Pop III galaxies. As a result, the Lyman- Werner background may never 
be able to build up to affect H2 molecules in minihalos. In combination, we suggest that 
cooling and hence star formation in minihalos can continue to take place largely unimpeded 
throughout the first reionization period, as long as gas is able to accumulate in them. 

II) We show that the intergalactic medium is likely to have been reionized twice, first 
at 2; = 15 — 16 and second at z = 6. Under a very conservative assumption for the star 
formation efficiency in minihalos with II2 cooling, the first reionization may be attributable 
largely to Pop III stars in large halos; in this case, processes promoting II2 formation in the 
cores of minihalos such as mentioned above in (I) do not really matter. In contrast, if star 
formation efficiency in minihalos is not more than a factor of ten less efficient than in large 
halos, then Pop III stars in minihalos may be largely responsible for the first reionization. 
In either case, it seems likely that Pop III stars in minihalos make a large contribution to 
the first phase of the metal enrichment of the intergalactic medium approximately up to a 
level where the transition from Pop III stars to Pop II stars occurs. 

This apparent inevitability of twice reionizations is reached, based on a joint constraint 
by two observational facts: 1) the universe is required to be reionized at 2; ~ 6 and 2) the 
density fiuctuation in the universe at z ~ 6 is well determined by the same small-scale power 
traced by the Lya forest observed and well measured at 5; ~ 3. As a result, the product of 
star formation efficiency and ionizing photon escape fraction from galaxies at high redshift 
is well constrained, dictating the fate of the cosmological reionization process. 

The prolonged reionization and reheating history of the IGM is more complicated than 
usually thought. Because both cooling and recombination time scales are much shorter 
than the Hubble time at the redshift in question, we devise a new, improved computational 
method to follow the evolution of the IGM in all phases, including HII regions, HI regions 
and partially ionized regions. The overall reionization process may be separated into four 
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joint stages. 

(1) Prom 2;~30tO2;~15 — 16, Pop III stars gradually heat up and ionize the IGM 
and the first reionization occurs at the end of this stage, when the mean temperature of the 
IGM reaches ~ lO^K. 

(2) The first stage is followed by a brief period with a redshift interval of order lS.z ~ 1, 
where the IGM stays completely ionized due to sustained ionizing photon emission by forming 
Pop III stars. During this period, the temperature of the IGM is maintained at ~ lO^K. 

(3) The transition from Pop III stars to Pop II stars sets in at the beginning of this 
stage [z ~ 13). The abruptly reduced (by a factor of ~ 10) ionizing photon emission rate 
causes hydrogen to rapidly recombine and the universe once again becomes opaque to Lya 
and Lyman continuum photons, marking the second cosmological recombination. Prom this 
time until z — ^ Compton cooling by the cosmic microwave background and photoheating by 
the stars sclf-regulate the Jeans mass and the star formation rate. The mean temperature 
of the IGM is maintained at ~ 10^ K from z ~ 13 to 2; = 6. Meanwhile, recombination 
and photoionization balance one another such that the IGM stays largely ionized during this 
stage with nnii/^^H > 0.6. Most of the star formation in this period occurs in large halos 
with dominant atomic line cooling. 

(4) At 2; = 6, the global star formation rate again surpasses the global recombina- 
tion rate, resulting in the second reionization of the universe. The second reionization is 
predominantly due to stars formed in halos where atomic line cooling is efficient. 

There is a wide range of interesting implications from this new reionization picture 
presented here. To highlight a few: 

• The Pop III stars enrich the IGM with a metallicity of ~ 10~^ Z0 at 2; ~ 13 — 15. 

• The magnetic field originating from massive stars could pollute the IGM with a large- 
scale coherent field (Z ~ lOOkpc comoving) of order ~ 10"^ G at 2; ~ 13 — 15. 

• The number density of Pop III massive black holes (10 — 300 M©) is comparable to 
that of globular clusters and should seed later structure formation. 

• Direct detection of Pop III hypernovae/supernovae are very difficult but a systematic 
search for gravitational lensing magnified Pop III hypernovae/supernovae targeted at massive 
clusters of galaxies may turn out to be fruitful; SIRTF may be able to detect them. SIRTF 
may also be able to detect some of the large Pop III galaxies. JWST should be able to detect 
both large Pop III galaxies and Pop III hypernovae. 

• The Thomson scattering optical depth is increased to 0.10 ± 0.03 (compared to 0.027 
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for the case of only one rapid reionization at z — 6), which will have significant implications 
on the polarization observations of the cosmic microwave background. Upcoming MAP 
results should be able to distinguish between these two scenarios. 

• Under the present scenario it would not have been a surprise that no Pop III stars 
have been found in the local universe. 

• The IGM, while opaque to Lya photons, is highly ionized {nm/nu ~ 0.1 — 0.2) up 
to redshift 2; ~ 7. This may reconcile the observation indicating the second reionization at 
2; ~ 6 with the detection of Lya galaxies at z > 6.5. 

• Finally, in contrast to the conventional reionization scenario, 21 cm lines prior to 2; = 6, 
which would otherwise be produced by minihalos in the absence of the first cosmological 
reionization, are predicted not to exist. This would provide a definitive test of the scenario. 

I thank Nick Gnedin, Zoltan Haiman and Martin Rees for stimulating discussion, Ren- 
nan Barkana, Ave Loeb and Jerry Ostriker for helpful comments and Rupert Croft for 
information on constraints on cosmological models placed by Lya forest observations. This 
research is supported in part by NSF grant AST-0206299. 
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A. Magnification of Pop III Sources by Gravitational Lensing 

We determine the fraction of area in the background source plane at redshift 2:5, that 
is magnified by a factor of > by foreground objects, primarily clusters of galaxies. The 
magnification regime in which we are interested is /x > 10. We will first derive a general 
relationship between the magnification (//) and source plane coordinate (x): fi oc 1/x for any 
singular density profile. Then, we will use the singular isothermal sphere thanks to its ease 
of analytic treatment, which has the same /i oc 1/x relation, to model clusters of galaxies 
to compute the gravitational lensing cross section at the high /i end. The computed lensing 
cross section is cross-checked by observations to ensure consistency. 

We express the magnification as a function of x, the source plane proper coordinate 
(Fig 13). It is well known that ~ 1/x in the small x, high /i limit, for both a point mass 
lens (Vietri & Ostriker 1983) and a singular isothermal sphere lens (SIS; Turner, Ostriker & 
Gott 1984). For simplicity we will consider only axisymmetric lenses. We will generalize this 
property of fi in the small ,t, high fi limit to any lens system as long as the lensing bending 
angle, 6, as a function of proper impact parameter, b, does not increase as fast as b. It is 
noted that 6* oc 6 for a uniform mass sheet, 9 = constant for a SIS, and 6* oc 1/6 for a point 
mass. From Figure 15 we can readily write down the lens equation as: 

« + = P-^: aDs^ODLs- (Al) 

The above equation relates the angular position of the image, a + (5, to the position of 
the source, /3 (both with respect to some reference direction) through the bending angle, 9. 
Ds, Dl and D^s are the angular diameter distance between the observer and the source, 
between the observer and the lens, and between the lens and the source, respectively, h is 
the proper impact parameter relative to the lens center. Expressing the bending angle 9 as 
9 = 9o{b/bo)'^ (m < 1), we can rewrite the above lens equation as: 

^o(V^'o)"^ = bDs/iD^DLs) - x/Dls. (A2) 

In the limit a; — >■ 0, the solution for the impact parameter is 6 = bo[boDs/{9oDLDLs)]^^^^''^^ 
(Einstein radius), and dx/db— (1 — m)Ds/DL. It is noted that the requirement m < 1 is to 
guarantee a solution, i.e., the existence of an Einstein radius. Then, it is straightforward to 
show that, at small x limit, the magnification is 

with the constant K, independent of x, being 

bpDs .i/im-i)Ds bo 
^9oDlDJ DLl-m- ^ ' 
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Fig. 17. — shows the gravitational lensing optics. O, L and S are the observer, the lens and 
the source, respectively. D^, Ds and D^s are angular diameter distance from observer to 
the lens, from observer to the source, and from the lens to the source, respectively, b and 
X are the impact distance and source plane coordinate, respectively. 6 is the light bending 
angle when passing by the lens at impact distance b. 
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We note that, for any (centrally) singular density profile, m is less than unity. It has been 
shown by Syer & White (1996) and Subramanian, Cen, & Ostriker (2000) that the density 
profile in the inner regions of dark matter halos, which are formed through hierarchical 
gravitational clustering/merging in the conventional Gaussian structure formation models, 
is p(r) oc where n is the linear power spectrum index at the relevant scales. For 

the scales of interest, we have n ~ —2 to —1 for cold dark matter models, giving p(r) oc r"^'^ 
to p(r) oc r~^^^. These inner slopes of halos were borne out in N-body simulations (Navarro, 
Prenk, & White 1997; Moore et al. 1999). Thus, we should expect // oc l/x at small x for 
halos formed in hierarchical structure formation models, including those of galaxy cluster 
size. Moreover, cooling and subsequent condensation of baryons in the centers of halos may 
further steepen the density profiles in the inner regions. 

Having deduced the universal /i oc 1/x relation at the small x end, we will now carry 
out further calculations by adopting the SIS model, because of its analytical simplicity and 
because of its astrophysical relevance as shown in many statistical studies (Gott & Gunn 
1974; Tyson 1983; Turner et al. 1984; Hinshaw & Krauss 1987; Narayan & White 1988; 
Wu 1989; Fukugita & Turner 1991; Mao 1991). This is further justified if /i is normalized 
at a somewhat lower value, where some data from current observations exist. We note 
that, for individual multiply imaged observed quasars, it is clear that one needs to include 
ellipticities of the lenses for realistic modeling. For example, one cannot produce quadruples 
with an axisymmetric lens (Narayan & Grossman 1989; Blandford et al. 1989). However, 
for our purpose of calculating the magnification cross section at moderate-to-high range 
/X ~ 10 — 100, the assumption of axisymmetry of lenses should be adequate. For a SIS the 
bending angle, ^o, conveniently independent of impact parameter, is (Turner et al. 1984) 

^0 = M-T, (A5) 
c 

where c is the speed of light and a\\ is the line of sight velocity dispersion of the SIS lens. 
With equation (A5) we can solve equation (A2) with the following solution at the x — > 
limit: b = and ^ = Inserting this solution into equation (A3) yields the small 

X limit of the x — /i relation for a SIS lens: 

x{n) = -i^fD^s. (A6) 
1^ c 

A SIS would subtend a sohd angle in the source plane AO, within which the luminosity of a 
source is magnified by > ii: Afi = ^^j^. Thus, the total probability of a random source at 
zs of being magnified by > /i can be obtained by adding up all the foreground lenses. This 
is done by integrating AQ/An over lenses of all masses and over all redshifts up to the source 
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redshift, zg-, resulting in a double integral: 

TTX^ dr 

Prandom^ J J -^^n{MA, z)4:7Tr^ — dzdM , (A7) 

where n(M^, z) is the mass function of halos at redshift z. To make the subsequent calcula- 
tion more analytically tractable, we assume 

n{MA,z) = {l + z)-^g{MA), (A8) 

where 

g{MA) = (4 X 10-VM,)(Ma/M,)-i [l + (M^/M,)"^] exp(-MA/M,) h^Mpc'^ (A9) 

is the differential cluster mass function at 2; = 0, taken from Bahcall & Cen (1993). Ma is the 
observed cluster mass within the Abell radius {ta = 1.5/i^^Mpc) and = 1.8 x 10^^ H^^Mq. 
The term (1 + is intended to approximately describe the evolution of the overall cluster 
mass function with redshift. The cluster mass, Ma, defined within the Abell radius, r^, may 
be related to a\\ for the SIS model by 

.H^, (MO) 

where G is the gravitational constant. Combining equations (A6) and (AlO) indicates that 
x'^ oc M|. It is instructive to examine the behavior of the following term, x'^g{MA), which 
contains all the dependence of the integrand of the integral in equation (A7) on Ma- Since 

x'^q^Ma) is constant at the low Ma end and x'^g{MA) oc (Ma/M^,) cxp(— M^/M*) at the high 
Ma end, the integral over M4 in equation (A7) is convergent on both ends and dominated by 
massive clusters near the exponential downturn at about M ~ 1 x IO^^Mq. Consequently, 
for our purpose, we need only to be accurate on the high mass end near the exponential 
downturn for g{MA) and its redshift distribution characterized by parameter w. Substituting 
equations (A6, AS, A9, AlO) into equation (A7) [and making use of the simple relations 
in an Qo = 1 universe: Dl = Rh{1 - l/v^TT^Z), Dg = RhC^ - l/V^ + zs), Dls = 
Rh{1/VT+^ - I/a/TT^) and ^ = |i?^^(l + z)-^f^] yields 



2B7t'G^R I 1 ^ 

random c^t*^ ^ 



2S 



2 

X 







x/TT^y V^TT^ \/rT^, 



1.0 X lO-^ETT^^G^iJ^ M2 / 1 r In 



-t)dt 



X 



1- / ) ( / - - j^ ) (1 + ^)-2/2-^d^. (All) 
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AU constants in equation (All) are in c.g.s. units, except Rh, which is in Mpc. Note that we 
have also inserted a constant B in equation (All), which serves to absorb the uncertainties 
due to other factors which either cannot be accurately treated here or are unknown, includ- 
ing background cosmology, deviations of density profiles from singular isothermal spheres, 
gravitational lensing due to other astronomical objects and uncertainties in the observed 
cluster mass function. The normalization of B will be set by comparing to observations at 
/I > 2. Now expanding all the constants and integrating equation (All) with respect to t 
(= Ma/M^), we obtain 

0.0365,, 1 ,_2,, , 

Prandom = 5" U " fr-, ^(^S, w), Al2 

11^ VI + Zs 

where I{zs,w) = /^^(l - -^)2(-^ - ^^^{1 + z)-^/^-^dz. The integral Iizs,w) can 
be done analytically (but the resultant expression is quite lengthy) and here we just give the 
final numbers for specific zs and w. For zs = 17.0, Prandom = (1-1 x 10^^,6.8 x 10~^,4.2 x 
10-^,2.1 X 10-^)5//i2, for w = (0.0,0.5,1.0,2.0), respectively. It shows that Prandom only 
weakly depends on w, because the integral is dominated by the moderate redshift range 
z~ 0.5- 1.0. 

It is justified to only consider the large splitting events since the lensing cross section 
is dominated by massive clusters which give rise to large splittings (> l"). For the sake of 
concreteness we will adopt w — 0.5, which is consistent with the relatively wild evolution 
of cluster density up to redshift about unity (e.g., Bahcall, Fan, & Cen 1997). As we will 
show below, the final results {Pdust) turn out to be extremely weakly dependent on w. 
Furthermore, we take B — 3.4 to normalize Prandom' 

2.0 X 10-^ 

i^randomxl^) — 5 ) l-^-'-'^j 

which gives Prandomil^ = 2) — 4.0 X 10~^. Since we are integrating the observed clusters to 
obtain the lensing probability and the same clusters are also responsible for observed gravi- 
tational lensing events, we can make a consistency check. For the Hewitt & Burbidge (1989) 
quasar catalog of 4250 quasars there are two confirmed multiple image lens systems with 
splitting > l", which corresponds to multiple lensing probability of 4.7 x 10~'^(= 2/4250), 
i.e., P = 4.7 X 10^'^ for /i > 2 (in the singular isothermal sphere case). Thus, our adopted 
normalization is consistent with the observed value of P = 4.7 x 10~^ for /x > 2. Considering 
possible selection biases (Turner 1980; Cen et al. 1994a; Kochanek 1995) and the counter- 
vaihng effect that the observed quasar sample is at significantly lower redshift than z — 17 
considered here, we argue that the adopted normalization is reasonable. From this we see 
that a random source at zs — 17 has a rather small probability of < 10~^ being very strongly 
lensed. 
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We will now examine the case where we do not sample the sky randomly, rather observe 
a set of selected regions centered on massive clusters. The cross section for magnification > 
in the source plane is 7rx'^{n), so the solid angle within which the sources will be magnified 
by > /i is 7ra;^(/x)/Dg. Integrating irx^lfi)/ over a pre-selectcd set of clusters with mass 
Ma > Miim in the redshift range z = zi to Z2 and dividing the integral by the number of 
clusters pre-selected gives the mean solid angle for source magnification > jj,. 

C B'g(MA)dM r 7^^ (1 + zy^Anr^^dz 

J^cl[Zl, Z2, Mlim) 

where Na{zi, ^2, Mu^) = g{MA)dM J^^{1 + z)-'"A7rr'^^dz is the total number of clus- 
ters selected (a normalization factor). Dividing the mean solid angle for source magnification 
> // by the field of view yields the mean probability of a source in such selected fields which 
will be magnified by > 

PciustW = —7^ , (A15) 

where Qpov is a telescope's field of view. As an example, let us take Mum = 1 x 10^^ Mq, 
zi = 0.0 and Z2 = 0.4, which yields A^^, (0.0, 0.4, 1 x 10^^) = (3020,2669,2362,1852), for 
w — (0.0,0.5,1.0,2.0), respectively, in an ^.q — 1, H — 50km/s/Mpc universe. Similar 
to equation (All), we can integrate equation (A14) for zs — 6.0, zi — 0.0, Z2 — 0.4 and 
Mum = 1.0 X lO^^Afo and the result is: Pdustilj) = (6.7 x 10"^ 6.7 x 10"^ 6.8 x 10"^ 6.9 x 
^^~'^^W^ 2b arcmin^ ^"^" adjustment parameter B' in equation (A14), which is 

different from B in equation (All) to reflect the fact that the uncertainties involved are 
somewhat different (for example, since the probability is normalized per cluster, the uncer- 
tainty in the amplitude of the observed cluster mass function cancels out), for clarity and 
simplicity we adopt B' = B = 3.3. The result indicates that Pdustil^) is nearly independent 
of unlike Prandom, which depends on w, giving 

6.7x10-^ Qpov y 

Pclust{^J) = — -2 — — — ■ (AI6) 

/i \25 arcmm / 

Comparing Pdust with Prandom, for the case with w = 0.5, for example, we see that one gains 
a factor of 34 in the cluster centered survey compared to a random survey, for a field of view 
of ^Fov = 25 arcmin^. The relative gain decreases with increasing size of the field of view, 
and the two probabilities become equal at ilpov =838 arcmin^, above which the calculation 
of Pdust is no longer valid due to multiple clusters in a single field. Note that, although 
the mean probability of a source within the field of view is inversely proportional to flpov, 
the total number of magnified sources per field of view is Pdust^FOV^src (where Hsrc is the 
surface density of source src, which is either gal or SN in this study) , which is independent 
of ^Fov- The physical reason for having a relative increase in the probability of magnified 
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sources is that the clusters of galaxies are rare targets and hence a random field of view has 
a small probability of intersecting a cluster. We note that Prandom is still valid even in the 
case of multiple clusters in a single field. However, the calculation of Prandom breaks down 
when multiple clusters are precisely aligned along the line of sight. But such cases are likely 
to be negligibly few. We further note that the spatial clustering of galaxy clusters would 
further boost the gain of Pdust over Prandom, but we will not consider this effect here. 



